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CONCLUSIONS 


1.  The  nature  of  natural  and  a pray  ralnfields  Is  such 
that  only  a  statistical  estimate  can  be  made  of  what  drops 
were  encountered  in  any  teat.  Hie  variability  of  the  slse 
distribution  of  drops  along  a  trajectory  can  be  appreciable. 

2.  Drop  sizes  around  0*5  mm  diameter  have  negligible 
water  content  and  can  be  Ignored  for  momentum  transfer. 

3.  For  impact  responding  fuzes,  the  larger  drops  from 
4  to  6  mm  are  most  effective,  and  the  test  facility  should 
characterize  the  number  of  such  drops  intercepted  by  the  fuse 
on  the  average,  and  their  standard  deviation. 

4.  The  question  of  how  high  a  rain  rate  can  be  used 
depends  on  the  time  needed  by  essential  fuze  components  to 
react  and  then  relax  sufficiently  to  an  impulse  encounter. 
Though  collision  time  with  a  raindrop  is  less  than  10  micro¬ 
seconds  for  weapon  firings,  the  response  of  the  fuze  depends 
on  its  configuration,  that  is,  the  stiffness  and  mass  values 
of  all  significant  discrete  elements  in  the  fuze.  However, 
there  is  some  reservation  in  trying  to  determine  these  times 
analytically,  as  the  analyst  may  omit  a  key  consideration, 
or  some  design  features  may  be  intractable  are lytlcally. 

Lach  fuze  needs  individual  consideration. 

5.  If  an  intense  spray  rainfield  is  chosen,  an  intensity 
of  at  least  20  in/hr  and  biased  with  larger  drops  should  be 
suitable.  There  should  be  sufficient  time  between  Impacts  to 
svold  the  question  of  multlpl  impacts  in  having  a  subsequent 
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drop  strike  before  the  fuze  has  effectively  relaxed  from  a 
previous  impact. 

6.  The  deluge  noizle  (80200)  used  at  Holloman  appears 
reasonable.  It  is  biased  toward  the  larger  size  drops. 

7.  The  O-C-O  (overhead  channel  orifice)  system  is  possible 
for  better  determination  of  the  drop  size  distribution.  However, 
horizontal  spacing  of  the  drops  is  needed  by  some  initial  motion 
in  the  channel  itself,  since  otherw.se  the  drops  would  be 
spaced  too  closely.  The  practicali  ,y  of  this  system  is 
questionable. 

8.  Statistical  uncertainty  is  practically  eliminated 

if  a  specific  drop  pattern  with  larger  drops  (principally  5  mm) 
can  be  positioned  on  a  planar  "web"  matrix,  normal  to  the  shell 
trajectory.  The  matrix  is  envisioned  as  a  thin  plastic  sheet 
with  discrete,  gelatinized  blobs  positioned  such  that  at  least 
one  of  them  viill  intercept  the  frontal  fuze  area  with  high 
probability.  (See  the  section  on  Rain  Web  for  details). 

The  drops  should  be  sealed  with  a  protective  coating.  Such 
sheets  can  be  positioned  so  at  least  20  milliseconds  of  travel 
time  exist  between  them.  This  deterministic  field,  while 
preserving  the  essentials  of  the  physics  of  the  encounter, 
will  provide  the  designer  with  an  idea  of  the  loads  to  be 
imparted  on  the  fuze. 

9.  Proximity  type  fuzes  (radio  frequency,  optical,  capa¬ 
citance)  require  further  consideration  of  the  effect  of  the 
test  site  on  the  fuze  electronic  response. 


10.  Aerodynamic  heating  would  soften  plastic  fuze 
frontal  structures.  Ordinarily  this  cannot  be  fully  simulated 
and  preheating  of  the  fuze  may  not  be  possible  because  of 
internal  explosive  leads. 


INTRODUCTION 


This  study  investigates  known  properties  of  rain  fields 
and  describes  technical  uncertainties  that  are  associated 
with  rain.  The  feasibility  of  establishing  a  military 
standard  must  address  such  technical  uncertainties  both  on 
the  part  of  the  rainfield  itself  and  the  response  of  the 
fuze.  The  literature  on  rain  phenomena  is  extensive  and 
information  that  pertains  to  such  a  standard  is  collected 
in  this  study. 

Topics  associated  with  rain  encounter  with  fuzes  are: 
rainfield  distributions,  for  both  natural  and  artificially 
produced  (spray)  rain;  reaction  time  of  fuzes;  accelerated 
testing  with  a  high  rainfall  rate;  the  practical  length  of 
a  test  facility;  and  the  method  of  obtaining  precise  rain¬ 
drop  distributions. 

This  study  addresses  principally  impact  initiated  fuzes, 
as  proximity  fuzes  (e.g .  RP,  optical,  capacitance)  are  also 
affected  by  the  proximity  to  the  ground  and  the  test  site 
structures. 

To  state  the  problem  succinctly:  What  test  arrangement 
is  most  feasible  from  an  economic  and  practical  aspect  that 
is  sufficiently  representative  of  the  physics  of  the 
encounter  of  an  impact  fuze  or  a  fuze  with  an  impaot  sensor 
with  rain?  Gan  the  test  facility  convey  a  predetermined 
loading  nattern  on  the  fuze  so  the  designer  C9n  anticipate 


loading  levels  on  the  fuze?  Can  a  sprinkler  type  system 
satisfy  these  goals  ?  What  are  the  least  number  of  fuzes 
that  can  be  fired  in  the  shortest  possible  rain  distance 
without  appreciable  statistical  uncertainty  about  which  and 
how  many  drops  intercepted  the  fuze  on  each  firing?  What 
maximum  rain  rate  should  be  allowed  for  accelerated 
testing? 

A  further  discussion  of  these  topics  is  delayed  until 
page  67  after  a  treatment  of  the  essential  elements  of  the 
rain  web  approach. 

A  previous  summary  study  on  rain  encounter  with  point 
detonating  fuzes  appeared  as  a  paper  in  the  Proceedings  of 
the  Fuze/Munitions  Environment  Characterization  Symposium 
(Picatinny  Arsenal,  Dover,  NJ,  November  1972)  and  is 
contained  as  appendix  B  in  this  report  since  it  serves  as 
a  basis  for  this  report.  (Reference  1) 


CONCEPTS  RELATING  TO  RAINFALL 


Measurement  of  Raindrop  Site  and  Number 

Various  drop  sizes  and  their'  number  that  are  instantane¬ 
ously  prerent  in  some  small  volume  are  measured.  After  many 
measurements  of  such  volumes,  the  average  number  of  drops 
according  to  same  drop  size  Intervals  is  reported.  This 
drop  size-number  distribution  is  normalized  to  a  represen¬ 
tative  volume  which  is  usually  one  cubic  meter.  (References  2,3) 
Questl ons : 

1.  /I hat  different  size  drops  are  in  the  sample  volume? 

2.  How  many  of  each  size  are  present  (some  size-diameter 

interval ' mus t  be  specified)? 

3.  What  is  the  spaclal  distribution  of  such  drops? 

4.  dhat  is  the  effect  of  wind  in  shifting  the  distribution? 

The  answer  to  question  3  might  be  a  Poisson  distribution, 

but  for  a  spray  rainfield  it  is  nighly  dependent  on  the  hard¬ 
ware  installation.  Question  4  may  be  treated  analytically  by 
wind  drag  factors. 

In  answer  to  questions  1  and  2:  Sizes  of  drops  vary 
from  a  fraction  of  a  millimeter  to  about  6  mm  diameter  (the 
very  large  drops  especiflly  may  not  be  very  spherical). 

Sizes  have  been  measured  in  various  ways,  from  sizes  of 
pellets  of  flour  formed,  by  photography,  and  to  expected  loser 
beam  measurement  of  shadows  produced  by  the  drops. 

Many  of  these  small  volume  measurements  are  needed  to 
obtain  a  representative  average  drop  distribution,  e.g, 
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if  the  sample  volume  were  as  large  as  100  cubic  inches,  about 
600  such  inspection  volumes  are  needed  to  obtain  one  cubic 
meter  of  representative  volume. 

Publications  that  report  raindrop  distributions  differ 
in  the  drop  diameter  interval  (ZD)  used  for  grouping  the 
range  of  drop  diameters.  Some  reports  use  0.1  mm,  0.2  mm, 
0.6mm  or  1  mm.  If  measured  in  0.5  mm  intervals,  the  number 
of  drops  of  mean  2.5  mm  diameter  actually  comes  from  drops 
in  a  range  from  2.25  to  2.75  mm  in  size. 

Reported  distributions  are  then  a  histogram  plot  of 
the  number  of  drops  of  a  particular  diameter  Dj  (within  a 

diameter  range  ZD)  as  shown  in  figure  1. 

Ni 


Figure  1,  Histogram  of  average  number  Ni 
of  drops  of  diameter  Di  for  a 
diameter  range  ZD  (here  0.5  ran) 
for  a  representative  volume  in  the  rainfield. 
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The  ordinate  is  plotted  on  a  log  scale  since  the  number 
of  various  sizes  changes  from  a  possible  several  thousand  small 
ones  to  a  fraction  of  drop  of  a  5  mm  diameter  in  a  cubic 
meter. 

The  total  number  of  drops,  N,  in  this  cubic  meter,  is 
not  too  meaningful.  It  .s  the  sum  of  the  drops  at  any  mean 
drop  size,  Di,  over  all  ,he  drop  diameter  increments, A D.  If 
the  drop  diameter  intervfl  is  large,  say  0.5  mm,  then  there 
are  fewer  intervals,  but  a  larger  number  in  an  interval. 

S  represents  the  number  of  discrete  average  drop  diameters  Dj: 

»  4  ,4  »i  (1) 

f  *  I 

Number  and  Size  in  Log-Linear  Distribution 

In  this  report,  in  order  that  the  reported  natural  and 
spray  ra infields  may  be  compared  In  a  convenient  manner,  all 
the  raindrop  distributions  were  resolved  to  the  same  AD  drop 
diameter  interval,  0.5  mm.  This  was  the  interval  used  In 
the  Holloman  report  of  1975  on  the  rocket  test  facility. 

Some  of  the  earlier  reported  distributions  used  0.1  nan 
intervals.  The  following  procedure  is  sufficiently  accurate 
in  converting,  for  example,  0.1  to  0.5  ram  intervals  by 
testing  the  water  content  of  the  transformed  distribution. 

This  approach  was  used  to  reduce  all  rain  distributions 
to  a  common  graphical  representation.  The  Synopsis  of  Back¬ 
ground  Material  for  MIL  STD  210B  (references  5  and  6) 
reports  distributions  on  a  1.0  mm  Interval,  and  only 
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a  crude  method  could  be  used  to  go  in  reverse,  to  a  smaller 
0,5  mm  interval. 

Water  Content  Comparison 


— 

1  .  1 

- —  - 

1  . 

)  - 

0  u 

NEW  0.5 
OLD  0.1 

Di-2flD 

Di-  A  D 

Di 

Di-*- A  D 

Dj  4  2  AD 

Old 

interval 

Pi 

D2 

d3 

°k 

P5 

1.8 

1.9 

2.0 

2.1 

2.2  mm 

sample 

Ml 

N2 

N3 

% 

H5 

Number 

drops 

The  mean  of  2.0  mm  diameter  will  be  used.  Change  the 
range  from  0.1  to  0.5  mm,  and  see  the  effect  on  the  liquid 
content,  which  goes  as  Nj.Di^.  The  liquid  content  in  tte  OLD 
system,  correct  to  terms  in  (Ad)^  ia; 


I  -  y  Pj  V;  =  0,  /v,  t-OjV*.  i-PjPs  +DV  K  +0*-^  (2) 

L  ou  o  . . 

U.-< D^f*  *'*-**’'•  +D>!  t(0'M0k,,C'HaClW 

i  -  fn  *-3  >, 

L<k-o  L  •»  , 

4.  +j  D*  (A0)  3^* 


(3) 


+ 

4  Id,’  *3  D?  +  3  p, 

4  [D53  4  3  P3^^0)4  3P,^of]/Vf 

The  llqui  1  content  of  the  eld  interval  system  is 
expressed  in  terms  of  the  liquid  content  of  the  new  with 
correction  terms: 
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L  *  cl  [ +  i  +/V3 
+  5D?ao[-  V*,  -Vi  *Pl 

4  J  DjCaof/V-  *l'*f'V*  +Wr3 


4-  Approximation  for  the 
liquid  content  for  the 
new  transformed  system, 

,  lnew. 


Correction  terns 
neglected  in  this  0.1 
to  0.5  mm  interval 
transfort*  tion. 


This  liquid  content  correction  is  small  and  depends  on 
what  drop  size  is  being  ionsidered  (D3  in  this  case;,  and 
the  specific  rain  distri  jution.  *‘or  example,  for  Marshall 
Islands  A,  "he  correctia  for  the  1  mm  drop  is  plus  of 
the  term,  LnkW.  *'or  the  2  ram  drop,  it  is  - 2 %  of  L^e^. 

dince  the  number  of  drops  changes  dramatically  from 
several  hundred  or  thousand  small  drops,  to  a  fraction  of  a 
drop  of  the  very  large  s  i.ze  in  a  cubic  meter,  several  orders 
of  magnitude  are  suitably  plotted  on  a  log  scale.  Natural 
rainfields  so  plotted  are  NOT  a  straight  line,  but  have  a 
maximum  around  2  mm  diameter.  For  natural  rain,  such  a  plot 
is  close  to  a  straight  line  from  about  2.5  to  5  mm  drops. 


Artificial  spray  fields  do  produce  a  large  number  of 
small  drops.  For  simplicity  of  presentation,  many  measured 
distributions,  as  for  the  Holloman  data,  are  reduced  by  a 
least  squares  fit  to  a  straight  line.  i'he  mathematical  form 
ofthe  equation  of  the  line  that  shows  this  distribution  is: 


1 


H 


Ni  =  No 

(-Di/Dr) 

e 

log  Ni  =  -1 

Dr 

[*] 

+  log  No 

(4) 


Ni  -  number  of  drops  in  a  cubic  meter  of  sice  Dj  (ram)  around 
a  range  diameter  & (mm). 

No  -  hypothetical  number  of  zero  diameter  drops  —  BOOO  AD, 
which  can  be  a  very  large  number. 

Dr  -  the  calculated  inverse  slope  of  the  line.  It  is 
sometimes  called  a  reference  diameter. 

There  is  reference  to  work  by  Marshall  and  Palmer  who 
related  this  diameter  to  the  rain  rate: 

Dr (mm)  =  0.244  R° * ^ ( mm/hr )  (5) 

The  drop  size  distribution  becomes: 

Ni  =  8000  AD  exp  {  -4.1  Di/R0,21)  (6) 


Liquid  Content  of  Rain 

Liquid  content  (grams  per  cubic  meter)  is  the  average 
liquid  volume  of  water  that  would  be  obtained  in  a  specified 
volume  (usually  cubic  meter)  cut  out  of  a  specific  ralnfield. 
Liquid  content  is  surprisingly  small,  even  for  heavy  rsins. 

^or  a  40  inches/hour  fall,  the  water  content  is  only  about 
0.004^  of  the  volume  of  the  rainfieJd. 

Definition  of  liquid  content  based  cn  the  measured  distribution 
L  (g/ra3)  -  WL  y  NiDi^  Z  IT  5*  NiDi3  (7) 

o  ‘V*  tooo  Z*  V" 
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with  the  dimensions:  Di  (mm) 

Ni  (drop3/cubic  meter) 
f>  -  one  gram/cc 

J-  number  of  discrete  drop  diameters 

Median  Drop  Diameter 

In  summing  the  liquid  content  of  all  the  drops,  that 
diameter  of  a  drop  at  which  of  the  liquid  content  lies 
below  this  diamet« r  is  called  the  median  drop  diameter. 

The  size  is  about  2  mm. 

Rainfall  Rate 

A  convenient  unit  to  use  here  is  inches/hour,  because 
by  a  numerical  coincidence,  the  rainfall  rate  and  liquid 
content  (g/m3)  are  about  identical  numerically  if  each  is 
expressed  In  the  units  given. 

Definition  of  rain  rate: 


J 

s 

R  (in/hr)  -  ViNlDi3  = 

>  ViNlDi3 
10»  c.— . 

■LU  ;*i 

(8) 

where  the  units  are:  vj^  (ft/sec)  terminal  velocity  of 

drop  of  diameter 

Di  ( mm ) 

Ni(  /m3) 

The  important  term  is  the  velocity  of  raindrops.  The 
higher  the  velocity,  the  higher  the  rain  rate  for  the  same 
drop  site-number  distribution. 

Acceleration  and  Terminal  Velocity  of  Drops 

One  reported  expression  for  the  measured  terul&al  velocity 
of  raindrops  is  a  quartic  polynomial,  where  the  velod  ty  Vj 
is  in  (m/sec),  and  the  diameter  Di  in  (mm):  (Reference  12) 

16 


(9) 


Vi  s  Cq  +  CiDi  C2Di2  +  C3Di3  +  C4D].1* 


The  coefficients  are: 

C0  *  -0.27128 

ci  =  5.22306 


02  "  -1.10757  CL  -  -0. 0014.6884 
C3  =  0.11115 


Table  1  expresses  formula  9  with  velocity  in  ft/sec: 

I  Table  1.  Raindrop  terminal  velocity 


D(mm) 

V(fps) 

0.5 

6.8 

1.0 

13.0 

1.5 

17.8 

2.0 

21.5 

2.5 

24.3 

3.0 

26.4 

3.5 

28.0 

4.0 

29.0 

4.5 

29.5 

5.0 

30.0 

5.5 

30.0 

It  is  profitable  at  this  point  to  go  one  step  further 

and  calculate  the  velocity  of  drops  as  they  accelerate  from 

rest.  There  are  two  reasons:  (1)  to  find  this  dynamic  fall 

under  gravity  from  a  liquid  stream;  (2)  to  find  the  horizontal 

induced  velocity  from  a  wind  gust  to  estimate  the  disturbed 

spacial  distribution.  The  velocity  of  spherical  liquid  drops 

starting  from  rest  under  gravity  is  calculated: 

Drag 


Motion 


i  © 

I 


Use  the  drag  force  expression: 


Drag  force  =  Cp  v  (K 

Cj)  s  experimental  drag  coefficient  (which  is 
a  function  of  Reynolds  number) 

»  air  densi  ty 

v  a  x  ■=  ins tantaneous  drop  velocity 

D  tdrop  diameter 

ms  mass  of  drop  a  fv.-^TT  d3/6 
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Equation  of  motion  of  drop  for  vertical  motion: 


m  x  =  mg-  CD  ±  V 

0 


Air-water  density  ratio  for  ordinary  conditions: 

0.00238/1.914-=  O.OCL2268 


The  acceleration  equation: 


•  • 

x  2.  g 


_  0.00092  CD 


(12) 


The  drag  coefficient  Cp  is  s  function  of  Reynolds  number  Re: 
,  x  Df*(A  _  x  D  _  x  D  (13) 


Re  = 


A 


0.00016 


where  the  Reynolds  number  itself  is  a  function  of: 


x-  Instantaneous  velocity  of  the  drop 
D~drop  diameter 
g.  air  density 
-  -dynamic  viscosity  of  air 

The  kinematic  viscosity  of  air  (^ki*  sA**^/  orc^^nary 

temperature  is  about  0,00016  ft^/sec.  Experimental  data  for 

spheres  in  air  is  from  Fluid  Mechanics  by  Binder,  p.  180.  The 
drag  coefficient  is  determined  from  the  current  Reynolds 


number:  Re 

CD 

Re 

CD 

1 

30 

200 

oTS  *) 

2 

18 

500 

0.58  l 

Range  of  interest 

9 

7 

1000 

0.49  f 

10 

4.5 

2000 

0.42  ) 

20 

3.0 

4-000 

0.40 

5o 

1.7 

100 

1.1 

These  equations 

were  solved 

for  drops 

from  1  to  6  mm 

diameter  to  find  their  time-fall' 

-velocity 

relation.  Figure  2 

plots  the  velocity  as 

a  function 

of  fd.  1  distance.  Ciroles 

represent  0.5 

second 

increments 

of  flight 

time.  Also  plotted 
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is  a  dotted  line  for  the  free  fall  of  any  object  with  no  air 
drag.  The  calculated  terminal  velocity  of  the  5  and  6  mm 
drops  exceed  the  experimental  value  of  30  fps.  Thi3  may  be 
due  to  treating:  these  drops  as  spherical,  when  in  nature 
they  are  not. 


Figure  2.  Spherical  drops  accelerating  from  rest. 
Spherical  Water  Drops  (Velocity  vs  Pall  height) 


Plgure  3,  Spberlcel  drops  at  0  and  8  fpa  Initial  velocity. 
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Figure  3  shows  superimposed  on  the  results  of  figure  2 
the  specific  case  when  drops  have  been  formed  from  a  stream 
which  has  already  reached  an  initial  velocity  of  8  fps. 

After  only  5  feet  of  fall,  drops  with  no  initial  velocity 
are  approaching  the  velocities  of  these  drops. 

These  plots  were  made  in  consideration  of  an  overhead- 
channel-orifice  (0C0)  system  where  the  channel  has  about  6 
inches  of  water,  with  an  exit  velocity  of  ^2gh  r 5.7  fps. 

In  experimentally  quiet  conditions,  about  one  foot  of  fall 
of  the  stream  is  needed  before  full  drop  formation  emerges. 

The  initial  velocity  is  set  at  about  8  ft/sec,  and  this  system 
is  discussed  In  a  later  section. 


Wind  Effect  on  Special  Distribution 


To  obtain  an  analytical  estimate  of  the  motion  of 
various  drops  for  a  horizontal  wind,  an  expression  similar 
to  equation  12  can  be  used.  Consider  drops  with  no  Initial 
horizontal  velocity.  A  constant  wind  velocity,  V,  is  applied. 
The  absolute  speed  of  the  drops  relative  to  the  ground  is  y. 
The  velocity  that  appears  in  the  acceleration  and  drag 


expressions  is  the  relative  velocity  (  V  -  y  ): 

i£j  *- 

\  h  =  0.00092  CD  (V  -  y)2/D  (Ik) 

~7T7^  Re  -  (V  -  y)D/»> 

The  consequent  horizontal  velocity  versus  horizontal  displace¬ 


ment  Is  obtained  for  two  cross  winds  of  1.55  and  I4..II  meters/ 


sec,  for  drops  from  1  to  6  mm. 
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Figures 4  and  g  show  these  trajectories.  The  total  length 
of  the  trajectories  Is  equivalent  to  the  time  these  drops 
should  fall  about  4  feet  vertically  in  air  under  gravity,  as 
expressed  in  equation  12.  Though  the  use  of  drag  coefficients 
for  spherical  objects  is  not  valid,  some  qualitative  idea  is 
obtained  of  the  sweeping  away  of  the  drops. 


SIDE  DISPLACEMENT  (CENTIMETERS) 

Figure  4.  Trajectory  of  stationary  drops  (1  to  6  mm  diameter) 
suddenly  exposed  to  cross  wind  of  3  knots. 


Figure  5.  Trajectory  of  stationary  drops  (1  to  6  bbb  diameter) 
auddenly  exposed  to  cross  wind  of  9  knots. 


By  table  1,  it  takes  an  appreciable  fall  for  larger 
drops  to  attain  terminal  velocity.  A  worse  condition  oocurs 
in  a  spray  ralnfield  if  the  nozzles  are  pointed  upward.  Their 
vertical  component  of  velocity  becomes  zero  at  the  top  of  the 
trajectory.  The  Holloman  ROCKET  track  nozzles  spray  downward, 
whereas  the  ARTILLERY  range  nozzles  spray  upward. 

If  a  specific  volume  is  cut  out  of  a  natural  ralnfield, 
and  the  drops  allowed  to  accumulate  at  the  bottom  of  this 
volume,  we  obtain  a  rain  gage  type  measurement  of  liquid 
content.  If  by  chance  v®  cut  out  the  same  volume  from  a 
3pray  ralnfield  and  find  that  the  liquid  content  Is  the 
same  (assume  drop  size  distribution  about  the  same),  the 
rain  rate  will  NOT  be  the  same. 


Spray 


Natural 


Drops  travel 
faster 


In  natural  rain,  the  drops 
are  travelling  at  a  higher  velocity 
than  In  a  spray  ralnfield.  The 
rain  rate  for  the  spray  ralnfield 
is  not  that  obtained  with  a  rain 
gage.  To  find  Its  equivalent  In 


nature,  the  drop  size  distribution  of  the  spray  ralnfield  is 


needed,  and  the  terminal  velocity  applied  to  the  various  drops. 
This  was  done  in  finding  the  rain  rate  of  such  rainfields 


with  experimental  terminal  velocities  in  equation  9. 
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PUBLISHED  NATURAL  AND  SPRAY  RAINFIELDS 


Introduction 


Table  2  lists  the  25  rainfields  Investigated.  Most  of 
them  are  from  reference  1.  Due  to  the  difficulty  in  ooraparing 
one  field  with  another  because  of  different  drop  diameter  range 
(4D)  used,  all  the  distributions  were  reduced  to  the  same  0.5  ran 
interval. 

Table  2.  Natural  and  Spray  Rainfields 

1  CANAL  ZONE  A 

2  CANAL  ZONE  B 

3  CANAL  ZONE  C 

4  CANAL  ZONE  0 

5  CANAL  ZONE  E 

6  STANDARD  RAINFALL  U.S. 

7  TROPICAL  RAIN 

8  MARSHALL  ISLANDS  A 

9  MARSHALL  ISLANDS  B 

10  MARSHALL  ISLANDS  C 

11  NORTH  CAROLINA  A 

12  NORTH  CAROLINA  B 

13  NORTH  CAROLINA  C 

14  MIAMI  CONTINUOUS  RAIN 

15  MIAMI  THUNDERSHOWER  A 

16  MIAMI  THUNOERSHOWER  B 

17  HOLLOMAN  ROCK  10-26-3 

18  HOLLOMAN  ROCK  10-26-5 

19  HOLLOMAN  ROCK  10-26-5SD 

20  HOLLOMAN  ROCK  10-26-6 

21  HOLLOMAN  ROCK  8070 

22  HOLLOMAN  ARTL  80200 

23  MILSTO-210B»MAX.ALT=0KM 

24  MRST0-21OB.  0.1*  MAX 

25  MILSTD-210B*  0.5*  MAX 


Distributions  1-5  are  the  heavier  rates  from  a  published 
series  for  the  Canal  Zone. 

Distribution  6  is  a  formulated  (theoretical)  rainfall  for 
the  U.S. 

Distributions  15-16  are  extremely  high  rates  recorded  In 
Miami. 
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Distributions  17,  18  and  20  are  Holloman  spray  rainfields 
in  the  rocket  test  facility,  measured  along  the  repeating 
horizontal  locations  from  10  to  26,  and  at  the  water  pressures 
of  3.5,  5,  and  6.5  psl.  The  designation  HOLLOMAN  ROCK  10-26-6 
refers  to  the  average  distribution  from  locations  10  to  26 
at  a  water  pressure  6,5  psi. 

Distribution  19  (HOLLOMAN  ROCK  10-26-5SD)  represents,  not 
a  real  distribution,  but  the  standard  deviation  of  measurements 
made  at  the  discrete  physical  locations  from  10  to  26.  It 
shows  the  appreciable  SPACIAL  variability  along  the  repeating 
8  foot  sections  of  the  rocket  traok. 

Distribution  21  (HOLLOMAN  ROCK  8070)  are  measurements 
around  1970  at  the  rocket  facility,  but  with  a  different 
nozzle  (8070)  which  produced  a  higher  rain  rate  of  about 
6  in/hr.  (Reference  2) 

Distribution  22  (HOLLOMAN  ARTL  80200)  are  measurements 
at  the  Holloman  artillery  (ballistic)  site  around  1969  with 
the  deluge  noszle  80200  with  rates  from  i*0  to  60  in/hr,  depend¬ 
ing  where  along  the  rainfield  measurements  were  made. 

Distributions  23-25  are  from  MIL  STD  210B,  and  will  be 
discussed  in  another  section.  It  was  difficult  to  fit  these 
in  as  they  are  reported  in  gross  1.0  mm  intervals.  To  achieve 
0.5  mm  increments,  the  number  of  drops  at  the  center  diameters, 
1,2,3...  was  considered  one  half  the  number  reported.  The 
number  of  drops  at  the  half  intervals,  1.5,  2.5...  are  calcu¬ 
lated  by  taking  one  fourth  of  the  drops  on  each  side  of  the 
associated  Integer.  The  calculation  of  rain  rate  and  liquid 
content  in  Table  4  is  done  with  these  0.5  mm  increment 
transformed  distributions,  and  with  equations  7  and  8. 


List  of  25  Fields,  Distributions,  and  Parameters 

Table  3  lists  the  average  number  of  drops  per  cubic 
meter  of  the  25  rainfields  considered,  at  each  diameter  from 
0.5  to  6  mm,  within  a  range  of  0.5  mm. 

For  example,  the  number  65  on  line  7  (Tropical  Rain) 
under  3.0  mm  diameter,  means  that  on  the  average,  65  drops 
should  be  expected  in  a  cubic  meter  for  this  rainfield,  for 
drop  sizes  from  2.75  to  3.25  mm  diameter. 

Note  the  spaces  for  Holloman  listings  17-20.  The  data 
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Table  3«  Drop  size  distributions  for  natural  and  spray  rainfielda 
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was  taken  from  a  1975  report.  The  Instrumentation  could  not 
read  drop  diameters  above  about  4.  mm.  Such  larger  drops  are 
expected,  and  Holloman  intended  to  update  the  instrumentation. 

Table  4  lists  typical  rainfield  parameters  calculated 
from  the  drop  distribution  in  a  computer  program.  The  total 
number  of  drops,  N,  in  the  last  column  is  not  too  meaningful 
because  of  the  large  number  of  very  small  drops,  which  have 
little  water  content  and  little  net  momentum  transfer. 


Table  4.  Calculated  Rainfield  Parameters 


NATURAL  and  spray  rainfields 

RATE 

(MM/HR) 

PATE 

(IN/MR) 

L IO 
(G/M3) 

MEDIAN 

(MM) 

N  DR0OS 
<  /M3) 

1 

CANAL  ZONE  A 

19. 

•  8 

.8 

1.8 

275. 

2 

CANAL  ZONE  0 

31. 

1.2 

1.3 

1.9 

377. 

3 

CANAL  ZONE  C 

51. 

2.0 

2.0 

2.1 

539. 

4 

CANAL  ZONE  0 

81. 

3.2 

3.2 

2.2 

884. 

5 

CANAL  ZONE  E 

117. 

4.6 

4.6 

2.1 

1464. 

6 

STANDARO  RAINFALL  U.S. 

SO. 

2.0 

2.1 

2.0 

2210. 

7 

TROPICAL  RAIN 

134. 

5.3 

5.1 

2.3 

1410. 

6 

MARSHALL  ISLANDS  A 

49. 

1.9 

2.2 

1.6 

1245. 

9 

MARSHALL  ISLANDS  0 

97. 

3.8 

4.2 

1.8 

2243. 

10 

MARSHALL  ISLANDS  C 

174. 

6.9 

7.0 

2.0 

2467. 

11 

NORTH  CAROLINA  A 

52. 

2.0 

2.2 

1.7 

936. 

1? 

NORTH  CAROLINA  $ 

115. 

4.5 

4.9 

1.8 

2456. 

13 

NORTH  CAROLINA  C 

191. 

7.5 

8.3 

1.7 

4448. 

14 

MIAMI  CONTINUOUS  RAIN 

688. 

27.1 

2  7.8 

1.9 

7614. 

15 

MIAMI  THUNDERSHOWER  A 

2344. 

92.3 

95.7 

1.9 

30664. 

16 

MIAMI  THUN0ERSH0WER  0 

2829. 

111. 4 

108.5 

2.2 

26914. 

17 

HOLLOMAN  ROCK  10-26-3 

56. 

2.2 

2.5 

1.7 

3540. 

18 

HOLLOMAN  ROCK  10-26-S 

57. 

2.3 

2.7 

1.5 

4389. 

19 

HOLLOMAN  ROCK  10-26-5SO 

29. 

1.2 

1.3 

1.7 

1626. 

20 

HOLLOMAN  ROCK  10-26-6 

57. 

2.2 

2.7 

1.5 

5125. 

21 

HOLLOMAN  ROCK  8070 

163. 

6.4 

7.5 

1.6 

6709. 

22 

HOLLOMAN  ARTL  80200 

1158. 

45.6 

41.9 

2.9 

18706. 

23 

MILST0-210B«MAX*ALT*0KM 

7780. 

306.3 

317.9 

2.1 

179683. 

24 

MILST0-210B.  0.1%  MAX 

495. 

19.5 

21.4 

1.7 

13672. 

25 

MILST0-210B.  0.5%  MAX 

100. 

3.9 

4.4 

1.6 

3027. 
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Plots  of  22  Fields ,  and  Comparisons 


Figures 6  to  10  plot  the  distributions  listed  in  table  2 


on  the  usual  log-linear  scale.  A  brief  summary  of  features: 


Figure  6.  General  trend  for  natural  fields:  As  rain  rate 
goes  up,  the  curve  shifts  up  rather  uniformly. 

Figure  7.  The  "standard”  rainfield  looks  like  a  spray  rain 
field  with  the  large  number  of  0.5mm  drops. 

The  tropical  rain  here  is  a  heavy  fall. 


Figure  8.  These  plots  of  Marshall  Islands  and  North  Carolina 
distributions  show  typical  upward  displacement  of 
the  curve  as  the  rain  rate  goes  up. 


Figure  9.  The  two  Miami  thundershowers  are  extreme  rainfalls. 
Note  the  appreciable  difference  for  the  larger 
drop  sizes. 


Figure  10.  Holloman  spray  fields: 


The  three  distributions  for  the  10-26  distance 
designation  are  about  the  same  in  rain  rate. 

As  the  water  pressure  increased,  the  number  of 
SMALLER  drops  increased,  as  the  number  of 
LARGER  drops  decreased. 

The  rocket  range  8070  nozzle,  and  artillery 
deluge  80200  nozzle  results  are  plotted. 


Comparison  plots: 


Figure  11.  Heavy  natural  rain  (TROPICAL)  is  quite  different 
from  the  1970  reported  Holloman  rocket  and 
ballistic  ranges  (8070  and  80200  nozzles). 


Figure  12.  TROPICAL  distribution  is  again  plotted. 

The  "standard"  field,  and  optimum  configuration 
of  Holloman  10-26-5  track  rather  well  with 
the  TROPICAL.  (Again,  no  data  on  large  size 
drops  from  Holloman). 

Figure  13.  'Aie  Holloman  deluge  80200  nozzle  Is  optimum 
for  creating  large  size  drops. 
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Rain  Distributions  in  Radar  Handbook 


In  the  Radar  Handbook  by  Skolnik  (reference  4)»  eight 
drop  distributions  for  precipitation  rates  from  0.25  to  150 
mm/hr  are  given  in  0.5  mm  diameter  increments  as  a  function 
of  percent  of  water  volume  from  each  of  these  drop  increments. 
The  book  cites  the  work  of  Laws  end  Parsons  reported  in  1914-3. 

To  see  how  the  distributions  differ,  compare  the  100 
mm/hr  ^3.94  in/hr)  rate  in  the  radar  handbook  with  Canal  Zone  E 
field  of  4.6  in/hr.  The  number  of  drops  per  cubic  meter  in 
the  handbook  is  obtained  by  using  the  rain  rate  equations: 

R  =  3.94  in/hr  s  7.2  TT^  vi  l  Pi  Di3/  106 
N  =  165  *' 

Ni=  N  Pi 

The  calculated  liquid  content  L  is  close  to  the  rain 
rate  in  inches  per  hour: 

L  =TT^Ni  Di3/  6000  5=  3.22  g/m3 
i 

Table  5  shows  the  comparison  when  the  percentage  of  liquid 
content  of  Canal  Zone  d  is  calculated. 


Table  5.  Comparison  of  Radar  Handbook  with 
Canal  Zone  E  Fain  Distribution 


Radar 

Handbook  1 

Canal  Zone 

E 

(3.94  in/hr) 

(4.6  in/hr) 

Size 

Percent  volume  Drops/m3 

Percent  volume 

Drops /m3 

(mm) 

Pi 

Ni 

Pi 

Ni 

0.5 

1.0 

1.7 

6. 

88 

1.0 

4.4 

7.6 

36. 

935 

1.5 

8.8 

14.7 

25. 

jfk 

2.0 

14. 

23 

17. 

250 

2.5 

17.1 

28 

9. 

134 

18.4 

30.5 

3.7 

55 

3.5 

15. 

25 

20 

4.0 

9. 

15 

0.42 

6.1 

4.5 

5.8 

9.6 

0.1 

1.44 

5.0 

_ 1±0 _ 

5. 

0.02 

0.32 

5.5 

2.8 

OI 

0.21 

6.0 

1.7 

- 

0.07 

6.5 

0.7 

1.1 

- 

i 


i 

l 

Though  the  rain  rates  are  similar,  the  distributions  are 
quite  different,  'i'he  radar  handbook  distributions  are  shifted 
toward  larger  size  drops. 

Rainfields  Reported  in  MIL  STD  210B 

MIL  STD  210B  (Climatic  Extremes  for  Military  Equipment, 
15  December  1973)  lists  environments  on  a  world  wide  basis 
(references  5,6).  Rain  distributions  are  listed  in  1  mm  dia¬ 
meter  intervals  from  an  analytical  expression  synthesized 
from  the  work  of  previous  authors.  Only  three  distributions 
are  given  In  table  6  for  sea  level.  These  respective  dis¬ 
tributions  are  listed  in  the  references  as  increasing  by  about 
30$  at  an  altitude  of  4  kilometers,  and  returning  to  about  the 
>  same  sea  level  value  at  around  8  kilometers  (4.6  miles). 

The  distributions  decay  to  about  half  the  8  kilometer  value 
rain  intensity  around  a  12  kilometer  altitude. 

Table  6.  Sea  level  and  8  kilometer  distributions  for 
Maximum,  0.1^  and  0.5^  probability  levels. 
(MIL  STD  21 OB) 


N«  s  Drops/me ter^ for  drop  - - > 


Rain 

type 

\  Rate 

Prec, 

diameter 

range 

given 

: 

(mm/hr) 

(in/hr) 

mass 

(g/m3 

0.5 

1.4 

1.5 

2.4 

2.5 

3.4 

3.5 

4.4 

4.5 

5.4 

5.5 

6.4 

MAX8 

1860. 

73. 

77. 

159,000 

30,010 

5,640  1064 

201 

38 

0.1$h 

188. 

7.4 

8.3 

11,800 

1,710 

247 

36 

5 

1 

o.5*c 

48. 

1.9 

2.2 

2,626 

342 

45 

6 

1 

1 

0 

Maximum  rate  ever  recorded  for  one  mimte. 

b  Rain  rate  should  not  exceed  this  0.1, <  Df  time  in 
wettest  part  of  the  world,  in  wettest  month. 

0  Rain  rate  should  not  exceed  this  0.5#  of  time  in 
wettest  part  of  the  world,  in  wettest  month. 


v 


I 
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*  VT-  'if  - 


A  correlation  was  found  between  the  rain  rate  R  and  the 


liquid  content  L,  and  expressed  as: 

L  (g/m3)  s  0.052  R  *  (mm/hr) 


(15) 


For  a  given  rain  rate,  the  pred.  pitation  mass  (liquid  con¬ 
tent)  is  calculated  from  equation  15.  This  is  the  precipitation 
mass  in  table  6.  To  go  one  step  further,  and  express  the  drop 
size  distribution,  first  a  median  drop  of  diameter  D0  was 
calculated  analytically,  again  from  the  rain  rate: 


Do  (mm)  1.48  R^*^ (mm/hr) 


(16) 


All  the  data  was  then  fit  to  a  log-exponential  function 
to  generate  the  original  drop  distribution,  Nj_,  for  a  drop 
diameter  interval  of  1  mm: 


where  Di =  1,2, 3, 4,5,6  mm.  Ni  are  numbers  reported  in  MIL  STD  210B. 

With  this  understanding,  the  numbers  for  drop  distributions 
in  MIL  STD  210B  are  obtained  roundabout,  through  analytical 
fit  of  original  data.  Care  has  to  be  used  in  using  these 
numbers  directly  because  of  the  coarseness  of  the  drop  diameter 
Interval  AD,  and  it  Is  not  expected  that  one  analytical  express¬ 
ion  should  fit  the  extreme  maximum  fall. 

Table  7  shows  the  liquid  content  for  the  three  fields 
calculated  with  equation  15>  which  is  a  valid  empirical 
expression,  and  with  the  general  liquid  content  expression, 
equation  5»  directly  from  a  distribution. 


Ni  =  389 


R1* 02 (mm/hr) 


-  3.67  Di/Do 


(17) 


34 


Table  7.  Liquid  content  (g/m3)  for  MIL  STD  210B 
Distributions  by  three  approaches 


L  -  0.052  R°*97 

T.  -  IT  + 

r  Ni°i 

Use  mid  point 

Di- 1,  2,  ... 

6000 

Use  lower  limit 

Di=  0.5,  1.5,  ... 

MAX 

77. 

lIt-6. 

340. 

0.1* 

8.3 

7.0 

18.5 

o.5* 

2.2 

1.4 

3.8 

Even  the  use  of  the  lower  limit  of  the  drop  sizes 
results  in  a  rain  rate  twice  what  it  should  be  for  the  MAXIMUM 
fall.  The  actual  liquid  content  for  the  0.1  and  0.5*  fields 
are  bracketed  by  the  two  choices  of  Di.  A  fi.  ner  reporting  of 
the  distributions  is  needed  to  work  with  them  directly. 

There  is  little  loss  in  accuracy  in  liquid  content  and 
rain  rate  In  going  to  a  coarser  diameter  Interval  ho.  if  the 
first  22  distribut. ons  In  this  report,  already  listed  for 
AD=  0.5  mm,  are  further  reduced  to  <1  Dal  mm  by  taking  the 
number  of  drops  at  the  integer  diameter,  and  adding  to  it  one 
half  the  number  at  both  sides  of  the  Integer,  the  resulting 
distributions  are  listed  in  table  8.  The  distributions  are 
distorted  at  the  smaller  drop  size,  but  this  Is  not  significant 
as  drops  around  one  millimeter  have  little  water  content. 

If  the  same  computer  program  is  used  to  reduce  these 
distributions  to  obtain  rain  rate  and  liquid  content,  the 
result  is  table  9.  The  overall  effect  is  to  artificially 
enhance  the  rate  by  about  9%.  In  the  program,  the  median 
diameter  Di-  1,  2,  ..  is  used  as  the  diameter  in  the  liquid 
content  and  rate  expressions  of  equations  7  and  8. 
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Table  8.  Distributions  reduced  to  1  ram  diameter  intervals 
for  comparison  with  MIL  STD  210B 


DROP  OIAMETER 

x 

z 

*** 

• 

o 

2.0 

3*0 

4.0 

5.0 

6.0 

i 

CANAL  ZONE  A 

118. 

138. 

16. 

i.l 

.0 

.01 

2 

CANAL  ZONE  B 

156. 

184. 

31. 

2.7 

.2 

.01 

3 

CANAL  ZONE  C 

230. 

240. 

53. 

7.2 

.7 

.07 

4 

CANAL  ZONE  0 

412. 

345. 

94. 

10.6 

1.0 

.12 

5 

CANAL  ZONE  E 

766. 

504. 

132. 

16.8 

1.1 

.18 

6 

STANDARD  RAINFALL  U.S. 

1204. 

209. 

48. 

7.6 

1.0 

.16 

7 

TROPICAL  RAIN 

685. 

498. 

145. 

24.0 

3.0 

.50 

0 

MARSHALL  ISLANDS  A 

737. 

407. 

27. 

1.0 

•  1 

0.00 

9 

MARSHALL  ISLANOS  B 

1278. 

616. 

91. 

3.2 

.4 

.11 

10 

MARSHALL  ISLANDS  C 

1306. 

849. 

194. 

17.4 

1.7 

.15 

u 

NORTH  CAROLINA  A 

533. 

350. 

42. 

3.5 

•  3 

0.00 

12 

NORTH  CAROLINA  B 

1600. 

650. 

92. 

12.5 

1.7 

.10 

13 

NORTH  CAROLINA  C 

2975. 

1185. 

144. 

16.0 

2.3 

.20 

14 

MIAMI  CONTINUOUS  RAIN 

2770. 

4040. 

743. 

50.0 

1.5 

0.00 

15 

miami  thundershower  a 

13620. 

13940. 

2225. 

180.0 

28.0 

2.00 

16 

MIAMI  THUNDERSHOWER  B 

11300. 

11780. 

2695. 

527.5 

80.0 

6.50 

17 

HOLLOMAN  ROCK  10-26-3 

2130. 

253. 

48. 

8.6 

0.0 

0.00 

10 

HOLLOMAN  ROCK  10-26-5 

2609. 

304* 

46  • 

5.9 

0.0 

0.00 

19 

HOLLOMAN  ROCK  10-26-5SD 

963. 

145. 

25. 

4.5 

0.0 

0.00 

20 

HOLLOMAN  ROCK  10-26-6 

2985. 

348. 

39. 

3.9 

0.0 

0.00 

21 

HOLLOMAN  ROCK  8070 

3875. 

1215. 

114. 

5.5 

.3 

0.00 

2? 

HOLLOMAN  ARTL  80200 

10125. 

2225. 

930. 

322.5 

90.0 

14.50 

If  the  rate  Is  calculated  for  MIL  STD  210B  distributions 
by  also  taking  the  integer  diameter  Di=l,  2,..  the  rain  rate 
is  highly  enhanced: 


Rate  (in/hr)a 

Rate  ( in/hr )D 

Maximum 

320. 

73. 

0.1^ 

16. 

7.4 

o.S% 

3.2 

1.9 

(a  Integer  diameter  in  eq.  8.  b  Reported  rate.) 

The  three  MIL  STD  210B  distributions  are  plotted  in 
figure  14  for  comparison  with  1)  tropical  rail,  2)  Holloman 
ballistic,  and  3)  Holloman  rocket  mean  field. 
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Table  9.  Calculated  rain  quantities  for  distributions 
reduced  to  1  ram  diameter  intervals. 


natural  and  spray  rainfields 

RATE 
( MM/HR) 

RATE 

IIN/HRI 

L 10 
I6/M31 

N  0R0PS 
<  /R3) 

1 

CANAL  ZONE  A 

22. 

•  9 

.9 

273. 

2 

CANAL  ZONE  B 

35. 

1.4 

1.4 

373. 

3 

CANAL  ZONE  C 

56. 

2.2 

2.2 

530. 

4 

CANAL  ZONE  D 

89. 

3.5 

3.4 

862. 

5 

CANAL  ZONE  E 

131. 

5.1 

5.0 

1420. 

6 

STANDARD  RAINFALL  U.S. 

60. 

2.4 

2.5 

1470. 

7 

TROPICAL  RAIN 

148. 

5.8 

5.5 

1355. 

0 

MARSHALL  ISLANOS  A 

58. 

2.3 

2.5 

1172. 

9 

MARSHALL  islands  b 

113. 

4.4 

4,7 

198b. 

10 

MARSHALL  ISLANDS  C 

196. 

7.7 

7.7 

2367. 

11 

NORTH  CAROLINA  A 

60. 

2.4 

2.5 

929. 

12 

NORTH  CAROLINA  8 

131. 

5.2 

5.4 

2356. 

13 

NORTH  CAROLINA  C 

22 2. 

8.7 

9.3 

4323. 

14 

MIAMI  CONTINUOUS  RAIN 

781. 

30.7 

30.6 

7604. 

15 

MIAMI  THUNDERSHOWER  A 

2651. 

104.4 

105.1 

29994. 

16 

MIAMI  THUNDERSHOWER  B 

3110. 

122.5 

117.0 

26389. 

17 

HOLLOMAN  ROCK  10-26-3 

70. 

2.7 

3.1 

2440. 

11 

HOLLOMAN  ROCK  10-26-5 

74. 

2.9 

3.5 

2964. 

19 

HOLLOMAN  ROCK  10-26-5SD 

36. 

1.4 

1.6 

1136. 

20 

HOLLOMAN  ROCK  10-26-6 

77. 

3.0 

3.7 

3375. 

21 

HOLLOMAN  POCK  8070 

202. 

7.9 

8.9 

5209. 

22 

HOLLOMAN  AR7L  80200 

1267. 

49.9 

46.1 

13707. 

The  object  of  this  exercise  is  to  show  that  although 
distributions  are  readily  plotted  for  comparison,  the  calcula¬ 
tion  of  rain  rate  and  liquid  content  for  good  accuracy  need 
a  AD  of  0.5  ram  or  smaller. 
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Figure  14.  Comparison  of  MIL  STD  210B  distributions  with 
three  others. 

a  MIL  STD  rain  rates  are  artificially  enhanced  if  the  integer 
diameter  is  used  in  the  rain  rate  equation  8.  The  three 
comparison  distributions'  rate  is  little  affected  by 
the  same  mathematics. 
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Cumulative  Water  Content  of  22  Ralnflelds 


Figure  1 5.  Cumulative  water  content  for  Canal  Zone  distributions. 

For  these  natural  rains,  as  the  rain  rate  increases,  the 
water  content  curve  moves  downward.  Although  the  drop 
distribution  curves  moves  rather  uniformly  upward,  the 
addition  of  larger  drops,  though  small  in  number, 
contribute  appreciably  to  delaying  the  curve  rise. 

Note  the  insignificant  water  content  for  smallest 
drop  sizes. 


Figure  16.  Cumulative  water  content  for  "Standard" 
and  tropical  distributions. 
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Figure  17.  Marshall  Islands  and  North  Carolina  again  show 
the  shift  downward  in  water  content  as  natural 
rain  rate  increases. 


Figure  18,  Cumulative  water  content  for  Miami  thundershowers 


DROP  DIRMETER  (MM) 


Figure  19.  Cumulative  water  content  for  the  Holloman 
rocket  and  artillery  (ballistic)  ranges. 

The  rates  for  the  thr>ee  curves  for  the  10-26  optimum 
locations  are  about  the  s^ime,  but  as  the  water  pressure 
Increases,  the  small  drop.!  are  increased  at  the  expense  of 
the  larger  ones.  Higher  pressure  for  these  noxzles  appears 
not  to  result  in  a  greater  efflux  of  water,  but  probably  a  more 
turbulent  outflow  to  tear  up  larger  drops. 

The  curve  for  the  Holloman  deluge  00200  is  unique. 

This  field  is  a  good  cand'date  for  assuring  encounters  with 
larger  drops  in  as  short  £  field  as  possible. 

Noticeable  for  the  ^clloman  ROCKET  facility  is  that 
already  at  1  mm  mean  diameter,  the  water  content  is  over 
25j6,  whereas  it  is  only  a  few  percent  for  heavy  natural  rains* 
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HOLLOMAN  AIR  FORCE  BASE  FACILITIES 
Rocket  Test  Track 

The  rail  track  is  50,000  feet  for  rocket  acceleration 
and  deceleration,  18,000  feet  of  this  length  is  available 
for  a  rain  field  (45  sections,  each  400  feet  long).  Each 
400  foot  section  has  50  nozzles.  Two  types  of  nozzles  were 
used  in  the  same  field  in  the  1975  report.  (References  7,8,9) 


Figure  20.  LATERAL  VIEW  of  nolloman  Rocket  track 


This  pattern  is  repeated  every  8  feet  as  shown  in  figure 
20.  Note  the  locations  10  to  26  where  rain  measurements  were 
made.  An  optimum  pattern  was  produced  along  this  line.  Nozzles 
are  positioned  about  6  and  7  feet  fromthe  track,  and  spray  down¬ 
ward,  about  60?  toward  the  track  as  shown  in  figure  21. 

Figure  21.  AXIAL  VIEW  of  holloman  rocket  track 
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In  1969,  a  single  type  8070  nozzle  was  used  on  this 
facility,  with  the  following  parameters: 

Nozzle:  VEE  JET  KLl+U  8070  (Standard^ 

31,000  optical  volumes  measured,  each  38  in^  (net  volume  19.2  m$) 
Rain  rate:  6.6  in/hr 
Liquid  content:  6.87  g/m3 

Median  diameter:  1.89  mm 
Nozzle  water  pressure:  9  pslg. 

Average  wind:  Cross  track  -  0.9  knots,  SD  1.5 
Along  track  -  1.3  knots,  SD  1.7 

Drop  si 78s  were  measured  to  about  ij.,5  mm  diameter. 

It  is  interesting  to  note  how  many  sample  volumes  are  needed 
before,  on  the  AVERAGE,  even  one  5  mm  drop  might  be  seen 
(size  4*75  to  5.25  mm).  The  inspection  volume  is: 

Volumes  38  in3s  0.000623  nr*  -  0.022  ft3 
With  the  8070  nozzle,  about  0.6  drops 

in  this  5  mm  range  might  be  present  in 

one  cubic  meter.  The  sample  volumes  needed 

before  on  the  AVERAGE,  one  such  size  Interval  drop  Is  met  is: 

1  Inspection  volume  m3  —  2700  sample  volumes 

O.OOO63  m3  0.6  drop 

This  represents  a  volume  of  60  cubic  feet  In  the  rainfield. 

If  a  fuze  wl  th  an  effective  fontai  diameter  of  1.5  cm  passed 

thraagh  the  entire  18,000  feet  of  rain,  thi3  is  equivalent 

to  only  about  half  this  volume.  Thus,  this  front  portion  of  the 

fuze,  in  two  flights  through  the  field,  would  on  the  AVERAGE 

meet  only  ONE  drop  in  the  interval,  4.75  to  5.25  mm.  This  is 

for  the  6  in/hr  rate,  8070  nozzlej 

Lower  (2.2  in/hr)  rate  nozzles  used  for  1975  report: 

Table  10  lists  the  rain  distribution  for  the  pptlmum  locations 
10  to  26  at  5  psi  pressure.  Note  the  spacial 
variation  within  this  8  feet. 
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Tabic  10,  Number  drops /m^  at  loeatlona  10  to  26 
at  5  pel  at  Holloman  roaket  traok 


LOC 

10 

a 

14 

16 

18 

20 

22 

24 

26 

ROW  AVG> 

SR 

43.1 

67.4 

38.7 

88.3 

101 

80.2 

21.2 

67.2 

.  61.6. 

62.1 

two 

1.92 

4.43 

1.98 

3.99 

4.37 

2.47 

1.21 

4.01 

2.70 

3.03 

MMD 

1.62 

1.43 

1.0? 

1,36 

1.36 

1.21 

0.82 

1.4) 

1.66 

1. 32 

N(O.S) 

1,120.0 

3,640.0 

3,820.0 

2, 990. 0 

3,020.0 

3,640.0 

3,230.0 

3,060.0 

1,330.0 

2.874.0 

M(l.0| 

498.0 

1,330.0 

763.0 

1,300.0 

1,420.0 

994.0 

573.0 

1,400.0 

714.0 

1,020.0 

Nd.S) 

217.0 

485,0 

181.0 

397.0 

490.0 

249. 0 

100.0 

474.0 

;94.0 

32I..C 

M(2. 0) 

104.0 

151.0 

58.7 

160.0 

222.0 

73.8 

28.3 

158.0 

139.0 

122.0 

N(J.S) 

34.9 

56.9 

27.5 

74.2 

83.3 

32.7 

11.3 

63.6 

47.6 

48.) 

W(J.O) 

13.3 

25.9 

7.79 

28.3 

26.0 

13.3 

1.38 

<3.7 

84.  T 

IS.) 

NO.  si 

3.23 

16.2 

1.86 

8.79 

12.1 

6,62 

0.46 

8.3) 

3.  S3 

T.02 

0) 

0.  44 

6.37 

1.88 

4.23 

1.68 

0.48 

0.94 

1.29 

Ml 

1.40 

SR  -  equivalent  rain  rate  (mm/hr) 

LWO  -  liquid  water  content  (g/m3) 

HMD  -  median  diameter  (non) 

N(D)  -  number  dropa/oublo  meter,  of  diameter  D  foraDsO,5wi 


figure  22  plots  the  number  of  drops/m^,  far  each 
drop  size,  as  a  function  of  the  nine  locations  along 
the  repeating  8  foot  spray  section. 

Figure  23  plots  the  average  number  of  all  the  drops, 
with  an  estimate  of  its  variation,  as  a  function  of  the 
drop  size. 
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Plgura  22  • 


Avaraga  Drop 
for  Drop  31 -a 
(Water  prea 


Artillery  (Ballistic)  Teat  Facility 


This  test  site  called  the  Hay  Draw  Test  site  is 
2000  feet  long,  and  shells  are  stopped  by  a  bunker.  Sprays 
are  positioned  on  opposite  sides  of  the  shell  trajectory. 
Figures  2k  and  25  illustrate  the  layout.  Data  on  this 
field  was  generated  with  the  deluge  noEzle,  80200. 


k-  50'  — J 

•  (Typical)  • 


-SPRAY  HEADS  ON  RISERS  (STAGGERED  *1'  APART) 
\JIARKER  TARGETS  (250'  APART) 

y  y 


INDIVIDUAL  RAIN  SECTIONS  (EACH  ROQ'  LONu 
RAINFIELD  (UP  TO  2,000'  LONG) 


?l£ura  24.  Layout  of  Sallisti°  Rain  Xast  Facility 


Update  to  ths  ballistic  range 

Nozzles  were  replaced  to  produce  a  rate  from  3  to  7 
in/hr  by  changing  pressure.  There  is  no  information  on 


47 


the  drop  size  distribution  with  these  nozzles.  Newly 
developed  instruments  are  expected  to  scan  shadows 


figure  25.  Ballistic  rainfleld  cross  section 


A  practical  consideration  is  shell  drop  under  gravity 
along  the  rainfield  length.  With  air  drag  neglected/ the 
drop  from  the  fcraJectory  top  is: 
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FUZE  ENCOUNTER  THROUGH  RAI'I 


► 


Mechanical  and  Electronic  Fusing 

A  partial  lists  of  Army  fuzes  is  given.  Several  have 
been  tested  in  rainfields  (references  10,  11), 


Artillery 

Mortar 

. Recoilless/Tank 

M557  PD 

M524A6  PD 

M91A2  BD 

m564  pd 

m?25  pd 

M503A3  PD 

M577  MT 

M567  pd 

M509A1  PIBD 

M732  PROX 

M734  MOP 

M530A1 

m534ai 

M539 

M578 

Rocket 

Missile 

Small  Caliber 

MI4-12E1 

M805 

M505  M579 

Ml+23 

M812 

M533  M594 

M427 

m8i£ 

M949  M714 

Ml+31 

m55o  mk  27 

Ml+33 

M438 

M551 

sal  velocities 

(fps):  Mortar  -  800 

Recoilless  ■ 

-  1200 

Artillery  - 
Tank  -  4-900 

1950 

Apart  from  the  frontal  erosion  effect,  which  is  partially 
treated  in  appendix  B,  fuzes  are  structurally  affected  in 
three  ways  in  high  velocity  rain  flight: 

1.  Single  impact.  Momentum  transfer  generally  depends  on 
the  cube  of  the  drop  diameter  and  linearly  with  the  shell 
qp  eed.  Larger  drops  are  effective  in  producing  high  stress 
levels  through  the  fuze,  depending  on  fuze  component  striffness, 

2.  Multiple  impact  (cumulative  effect)  occurs  as  in  a  pro¬ 
gressive  crushing  or  subsequent  effective  impacts  before  the 
fuze  has  sufficiently  relaxed.  This  occurs  usually  after 
some  critical  drop  diameter,  and  if  the  drops  are  too  closely 
spaced  as  might  occur  in  an  accelerated  test, 

1+9 


3.  Resonance  effect  in  small  supported  structures,  and 
microphonics  in  proximity  and  electronic  fuzes.  If  there 
is  some  periodicity  in  the  rainfield  (or  if  a  driving  oscilla¬ 
tion  is  from  some  aerodynamic  effect  a  lone ) ,  the  fuze  compon¬ 
ents  will  react  at  this  frequency,  with  a  larger  amplitude  If 
the  periodicity  is  near  the  Resonant  frequency  of  the  fuze 
components. 

Fuzes  are  too  complex  generally  in  their  reaction  to 
permit  a  detailed  analytical  look  by  wsy  of  structural  response 
to  all  possible  reacticr  s  to  impulsive  frontal  Impact. 

Specific  Case  with  Tropical  Rainfield 


As  an  example,  con.-d  der  a  fuze  with  an  effective  frontal 
diameter  of  1.5  centimeters,  and  that  drops  hitting  this  area 
of  1.767  cm2  will  affect  the  fuze.  Momentum  transfer  to  the 


lateral  surface  of  the  ogive  Is  degraded  and  Is  not  considered. 
The  shell  travels  in  a  100  meter  rainfield,  with  a  velocity  of 


1000  meters/sec. 
1.5  cm 


A*  area  2  1,767  cm2 

Ls  rainfield  length  2  100  m  «  328  ft 

v*  shell  velocity  *  1000  m/sec  r3280  fps 


ACTUAL  Consider  th? ee  AVERAGE  parameters  of  this  encounter: 
SIZE 


1. 


the  distance  In  meters  the  shell  moves,  ON  THE 
AVERAGE,  before  it  Impacts  one  such  raindrop: 


1 

A  Ni 


2.  the  number  of  Impacts  with  these  drops  made  ON  THE 

AVERAGE  In  shell  travel  through  100  meters: 


Pi  =  L  A  Ni 
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3.  ATi  the  time  (seconds)  between  Impacts  ON  HIE  AVERAGE 
of  such  size  drops  on  the  fuze  frontal  area: 

A  Tj  =  At-  -  - i - 

v  **  A  v  Hi 

For  this  example,  the  TROPICAL  field  of  132  mm/hr  (5.3 
in/hr)  is  used.  Ni  is  the  number  of  drops  of  diameter  Di 
(40  =  0.5  mm)  per  cubic  meter.  Table  11  lists  the  AVERAGE 
encounter  values. 


Table  11.  AVERAGE  encounter  values  in  TROPICAL  FIELD  (5.3  in/hr) 


Di 

(mm) 

Ni 

(  /m3) 

Ai  z  a  N; 

Distance  (meters) 
before  one  drop 
encountered 

-  L  A  N; 

Number  impacts 
for  100  meter 
travel 

A  rJ~  N; 

Time  (seconds) 
between  drop 
impacts 

1.5 

360 

16 

6 

0.016 

2.0 

250 

23 

4 

0.023 

2.5 

135 

k2 

2.4 

0.042 

65 

87 

1.1 

0.087 

25 

226 

0.44 

4.0 

10 

566 

0.18 

KBzSHB 

4.5 

3 

1890 

0.053 

5.0 

1.25 

(4-530 

0.022 

1 

5.5 

■Bi 

0.009 

11.0 

6.  0 

Kb 

I 

0.004 

23.0 

Observations  on  this  simple  example: 

1.  Pi,  the  number  of  Impacts,  deperd  s  on  the  field 
length.  Field  length  increase  Is  done  at  the  expense  of 
shell  drop  and  maintaining  uniformity  of  rain  over  the  tra¬ 
jectory.  For  a  1000  foot  field,  multiply  these  impacts  by  3. 

2.  In  all  three  parameters,/;,  Pi,ATj.,  the  actual  rain 
distribution,  Ni,  Is  present.  For  fuse  response,  It  is  not 
sufficient  to  Just  specify  rain  rate,  but  the  rain  DISTRIBUTION, 
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to  ascertain  the  average  number  of  large  size  drops  present, 

3.  The  important  factor  is  AT.  For  a  given  fuze  and 
velocity,  how  intense  can  the  field  be  made?  Statistics  can 
readily  be  used,  but  for  now  consider  average  values  of  AT, 

Figure  26  shows  this  TROPICAL  distribution  increased  by 
a  factor  of  I4.  and  8,  The  average  time  AT  becomes  divided  by 
these  factors.  For  a  factor  of  U  (about  20  in/hr),  the  average 


Figure  26.  increase  of  TROPICAL  distribution  by  fsetors 
of  4  and  8. 

time  between  collisions  for  drops  greater  than  3  mm  is  not 
less  than  20  milliseconds.  This  may  be  a  sufficient  time 
for  fuze  response  relaxation. 
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Again,  with  this  TROPICAL  field  increased  by  an  intensity 
of  4  to  about  20  in/hr,  for  a  4.5  mm  drop  (0.5  mm  interval), 
with  the  average  time  between  collisions  for  these  drops 
now  about  470  milliseconds,  about  1500  feet  of  such  a  20  in/hr 
rate  would  have  to  be  traversed  for  the  1,5  cm  diameter  fuze 
tip  to  engage  ONE  such  drop  (on  the  AVERAGE). 

The  point  is:  Are  a  few  larger  or  many  smaller  (or  both) 
drops  more  detrimental  to  the  fuze?  If  the  larger  ones,  then 
increase  the  rain  rate  as  high  as  possible,  and  even  bias  it 
with  one  with  many  large  drops. 


RAIN  WEB  TOR  FUZE  RAIN  TESTING 


Introduction 

Some  evident  handicaps  of  rain  testing  for  PD  fuses  are 
1)  wind  gusts,  2)  elaborate  physical  layout,  3)  uncertainty 
of  type  of  drops  intercepted  in  any  firing  with  the  need  of 
many  firings  to  lend  credibility  to  the  fuze  response.  Effect 
of  aerodynamic  heating  must  necessarily  be  missing.  These 
capricious  aspects  suggest  the  use  of  a  rain  web  approach  as  a 
compromising  solution. 

The  UK  has  experimented  with  plastic  spheres  in  excess 
of  5  mm  diameter  and  with  captive  rain  drops  suspended  in  a 
crosshatch  web  system.  Gun  firings  have  been  done. 

Outlined  here  is  an  example  of  a  "rain  web"  which  is  meant 
to  obviate  these  handicaps.  While  the  essential  feature  of  rain 
encounter  is  preserved,  this  is  accomplished  in  an  arrangement 
that  is  convenient,  easy  to  Implement,  and  has  a  predictable 
distribution  of  drops  intercepted  for  any  firing. 

Since  larger  diameter  drops  from  ^  to  6  mm  are  important  for 
impact,  a  planar  array  of  5  mm  diameter  drops  will  be  usedfor 
this  study  as  an  illustration. 

Fabrication  of  the  Drops 

In  some  manner  the  drops  are  held  in  a  stationary  pattern. 
They  can  be  gelatinous  blobs,  with  additives  to  prevent  freezing, 
and  contained  between  thin  sheets  of  plastic  perhaps  sealed  by 
heat.  If  possible,  to  keep  weight  down,  the  drops  might  be 
sealed  individually  or  sprayed.  Depending  on  the  substrate 
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strength,  holes  night  be  initiated  between  drops.  Details 
depend  on  material  properties. 

If  all  the  drops  are  the  same  size,  one  repeating  pattern 
is  established  with  13  drops,  where  the  horizontal  center  to 
center  spacing  is  a  distance  L,  so  there  are  two  drops  within 

a  cell  in  area.  (Actual  closest 
center  to  center  distance  is  L  sin  45°). 

A  water  drop  of  5  mm  diameter  has 
a  volume:  TTd3/6  s  0.0655  cm3. 

Consider  for  this  example  a  fuze  with 
an  effective  frontal  diameter.  Dp  of 
1.5  cm.  If  a  plastic  sheet  supporting  the  web  is  0.002  inch 
thick,  (assume  a  density  of  water  for  the  plastic),  the  volume 
of  sheet  intercepted  by  the  fuze  front  is: 

7T  Dp2(0.002)(2.54)  =  0.009  cm3 
— - u - 

In  this  example,  this  plastic  portion  intercepted  by  the 
fuze  tip  represents  14 £  of  the  volume  of  one  of  the  drops 
and  stresses  the  need  of  keeping  the  supporting  web  weight 
low,  as  this  delivers  a  momentum  over  the  entire  frontal  area, 
not  particularly  representative  of  rain  encounter. 

If  the  volume  in  •  5  mm  drop  were  transferred  to  an  ideal 
cylinder  of  base  either  4  or  5  mm,  the  height  would  be: 

tj*  3  mm 

mm 


4—  L— * 
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Though  refinements  to  the  computer  program  can  include  the 
effect  of  the  weight  of  the  sheet,  and  of  non-uniform  drop 
cylinder  forms,  this  simulation  will  consider  5  mm  drops  as 
5  mm  diameter  cylinders,  4.1  mm  high,  with  a  weightless 
support. 

It  is  well  to  estimate  the  water  weight  (pot  including  the 
weight  of  the  supporting  web)  of  a  3x3  foot  web  structure  for 

5  mm  diameter  drops  with  a  spacing  L=15  mm  (closest  center-to- 

0 

center  distance  would  be  L  sin  45=10.6  ram). 

Weight(lb).  W2  (?.54)2tT03  2  (12)2 

~  60  L 2  (453.6) 

0.2145  (W  ft)2  (D  mm)3 

»  - UtL -  ' -  —  1.07  lb  of  water. 

(L  mm) c 


Frequency  Distribution  of  the  Encounter 

In  the  example,  all  possible  orientations  of  the  fuze 
with  the  web  are  considered  if  the  center  of  the  fuze  frontal 
area  goes  through  a  random  selection  of  points  within  the 
rectangle  of  area  L/4  *  L/2  defined  in  Figure  27: 


Outer  circumfer-  CQ 
ence  of  fuze 
frontal  area 
with  center 
located  at 
lower  left  end 
grid  point  of 
the  rectangle 
engagement  grid.  ~ 


«t.j  i 


; /gyi  0 

©  © 

©  © 


Figure  27. 

Basic  13  drop 
pattern  for  Rain 
Web  with  engage¬ 
ment  rectangle. 
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The  basic  supposition  is  that  the  approaching  center  of 
the  fuze  tip  can  pass  through  any  point  on  the  rain  web. 

The  grid  pattern  for  both  the  rain  web  and  the  approaching 
fuze  was  done  in  the  following  manner: 

1.  A  repeating  13  drop  pattern  was  represented  by  a  I6lx  161 
point  grid,  or  25 921  grid  points.  Zeros  at  the  points  represented 
locations  of  no  drops,  and  l's  where  the  drop  areas  were  present. 
An  x-y  coordinate  system  has  its  origin  at  the  center  of  drop 
number  5. 


10  cp 

©.-I-..  © 

ife. 


*/ 


I 


l+2. 


The  distances  involved  are: 

Da  4.5  mm  (representing  5  mn 
diameter  drop). 

Dp=  15  mm  (fuze  diameter). 

L  varies  from  11  to  23  ran. 

2.  Locations  of  the 
boundary  and  inner  grid  points 
of  che  13  circles  were 
determined  by  the  inequality: 


where 

the 

centers 

of  the 

13 

circles 

are  located 

at  A 

A 

B 

A 

B 

A 

B 
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-L 
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L 

0 

10 
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L/2 

2 

0 
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7  ■ 
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-L 

11 

L/2 
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3 

L 

L 

8 

0 

-L 

12 

-L/2 

-L/2 

4 

-L 

0 

9 

L 

-L 

13 

L/2 

-L/2 

5 

0 

0 

In  the  simulation,  each  of  these  13  eouations  was  in¬ 


terrogated  by  sweeping  through  the  25921  combinations  of  x  snd  y 
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values.  When  this  inequality  was  satisfied,  th9t  particular 
grid  point  (x,y)  represented  a  region  inside  the  drop. 

3.  The  number  of  distinct  encounters  or  collision  orientations 


of  the  fuze  with  the  web  is  the  number  of  grid  points  inside  the 
collision  rectangle,  which  is  essentially  the  area  ratio  of 


this  rectangle  to  that  of  the  entire  square: 


Number  of  distinct 
grid  encounters 


'fi/i'lQ./v)  I  . 
[jf/.+j)!'. 


558  for  L  a 11  mm. 
612  for  L*  15  mm. 
646  for  L  -  19  mm. 
67 2  for  L  =  23  mm. 


4.  The  fuze  tip  is  presented  as  a  circle  with  the  center 
located  at  the  coordinates  W  and  V.  As  the  fuze  engages  the 

web,  the  center  of  the  fuze  tip 
located  at  V,W  can  be  at  any  of 
the  648  grid  points  inthe  rec¬ 
tangular  array  Cfor  L=19  mm  spacing). 
The  location  of  the  fuze  area  then 
is  determined  by  interrogating 
the  entire  161x161  grid  points  of  x  and  y  for  specific  cases 
of  V  and  W.  Those  values  of  x  and  y  which  satisfy: 


are  grid  points  where  the  fuze  tip  has  casts  its  profile  on 
the  rain  web. 

5.  Counting  raindrop-fuze  tip  overlaps: 

Consider  the  L=19  mm  spacing.  For  the  possible  648  encounters, 
the  background  161x161  gi id  was  interrogated  for  both  the  drops 
and  the  fuze  tip.  Where  an  overlap  occurred  with  the  fuze  tip 
and  any  portion  of  a  drop,  a  1  was  assigned  to  this  grid  point. 
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For  each  encounter,  the  drop-fuze  overlap  was  counted  and 
divided  by  the  number  of  l's  which  constituted  one  of  the  drops 
in  the  grid,  to  obtain  the  effective  number  of  5  nun  drops 
encountered.  Fractions  of  drops  were  summed  to  obtain  the 
effective  number  of  drops. 

For  example,  for  0ne  orientation,  520  overlapping  l's 
were  counted  between  the  fuze  and  the  drops  it  overshadowed. 

In  this  case,  3l*l  grid  points  represented  a  cylindrical  drop. 
This  constituted  hitting  520/341-  1.525  drops  of  diameter  5  mm. 

To  find  the  probability  of  encounter  with  1.525  drops, 
the  total  number  of  times  that  520  overlapping  l's  occurred 
for  the  6I4.8  possible  orientations  was  calculated.  The  final 
sum  Indicated  this  occurred  7  times,  and  the  probability 
became:  PDl  =  1 .  c;2e;  =  *  0.0108 

h'our  examples  are  given  with  uniformly  distributed  5  mm 
drops  (here  represented  as  i| .5  mm  diameter  cylinders),  as  the 
drop  center  separation  L,  as  defined  in  Figure  27,  varied 
as  11,  15,  19,  23  mm.  The  fuze  diameter  remained  at  15  ram. 

Figure  28  shows  this  array  with  L -  11  mm  (a  rather  close 
drop  spacing),  and  the  contour  of  the  fuze  diameter  centered 
at  the  lower  left  end  of  the  L/kxL/2  encounter  rectangle. 

After  sweeping  the  entire  rectangle,  the  program  g en era  ted 
the  freauency  distribution  for  interce -sting  the  number  of  5  mm 
drops  on  any  fir'ng.  This  is  Figure  29,  where  the  frequency 
distribution  is  normalized  to  one.  (Without  normalization,  the 
area  under  this  curve  is  one).  In  any  firing  with  L = 11  mm, 
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there  53  a  high  probability  of  intercepting  2.5  drops,  with 
none  below  this  number.  The  distribution  goes  to  zero  just 
before  I4.  drops.  The  probability  of  interception  for  a  number 
range  is  estimated  from  this  plot  by  taking  some  width  along 
the  horizontal  axis,  and  noting  the  area  above  this  width 
relative  to  the  area  under  the  entire  curve. 

The  freouency  distribution  was  somewhat  grainy  from  the 
181  x  181  grid,  so  a  least  sauare  polynomial  matrix  smoothing 
was  used  with  a  linear  fit  over  several  grid  points  to  obtain 
the  representative  distribution  of  Figure  29.  An  excessive  num¬ 
ber  of  points  was  avoided  to  prevent  distorting  the  distribution. 


Figure  30  is  the  L=  15  mm  separation  of  drops,  with 
Figure  31  as  the  resulting  frequency  distribution  of  intercep¬ 
tion.  The  range  is  now  from  one  to  two  drops,  in  a  rather 
unioue  distribution. 

Figure  32  is  for  L  -  19  mm  separation,  and  Figure  33  shows 
a  rather  convenient  distribution  peaked  at  one  drop,  decreasing 
to  zero,  with  the  other  side  going  to  1.35  drops. 

The  actual  physical  si7e  of  this 
seemingly  optimum  array  is  given 
at  the  right. 

Figure  3I1  is  the  distribution 
with  L  increased  to  23  mm.  The 
peak  occurs  at  one  effective 

drop  intercepted,  but  now  the  probability  of  collision  with 
less  than  one  and  even  no  partial  collision  at  all  increases. 
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Figure  28.  Web-fuze  configuration  for  L  -  11  mm. 

Dashed  circle  represents  circumference  of  l5mm  diameter  fuze. 
Dotted  rectangle  represents  encounter  points  with  fuze  carter. 
(This  drop  configuration  is  very  dense). 


S 


Figure  29.  Calculated  drop  interception  distribution  (L—  11  mm). 

(Portions  of  drops  are  summed  to  the  effective  number 
intercepted ) . 
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Figure  30.  Web-fuze  configuration  for  ram. 

Circles:  k»S  ram  cylinders  (5mm  drops) 

Dashed  circle:  15  ram  fuze  front 

Dotted  rectangle:  Sample  space  of  collision  orientations 


EFFECTIVE  NUMBER  OF  5  MM  DROPS  INTERCEPTED 


Figure  31 •  Calculated  drop  interception  distribution  (b =  15  ram) 


NORMAL  I 
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^  g  i.uu  i  «uu  j»uu  s 

EFFECTIVE  NUMBER  OF  S  MM  OROPS  INTERCEPTED 


Figure  33.  Calculated  drop  interoepfcion  distribution  (L»  19  an) 


63 


Figure  32.  web-fuze  configuration  for  L e 19  mm. 

Circles:  mm  cylinders  (5mm  drops) 

Dashed  alrcle:  15mm  fuze  front 


S 


effective:  number  of  s  mm  drops  INTERCEPTED 


Figure  33.  Calculated  drop  Interception  distribution  (Ls  19  mm) 
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Figure  34.  Calculated  drop  intercept  ion  distribution  (L- 2  3  ran) 

(Increased  drop  spacing  has  increased  the  probability 
of  no  partial  collisions.) 


Fractional  drops  intercepted  and  Water  Content 


•‘■he  sane  res'ilts  should  be  obtained  (and  were)  for  a  simple 

9  drop  souare  array,  where  the  collision  encrunter  square  ia 

2 

now  (L  sin  45°/2)  on  a  side:  Area  -s  (  L  sin  45\2  _  L_ 

x  2  /  “  8 


,L  ain  45 


©  0  .0 


i 

L 
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The  analysis  was  extended  one 
step  further  to  find  the  distri¬ 
bution  of  fragments  of  drops  and 
how  much  of  the  drop  encounter 
comes  from  impacting  an  integral 
drop,  and  how  much  of  the  encounter 
comes  from  adding  sums  of  several 
drops  impacted  to  sum  to  one  drop. 


The  drops  were  numbered  and  for  any  collision  orientation,  the 
l»s  for  each  drop  were  relegated  to  0  histogram  of  drop  portions 
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in  0.025  volume  segments  of  a  drop 


Only  9  drops  are  needed  for  the  range  of  parameters  used 

A 

in  this  example.  The  encounter  square  now  had  (JU) 

s  35x35-1188  collision  orientations  for  the  L=19  ram  spacing. 

Figure  35  plots  the  normalized  volume  of  water  expected 
for  any  encounter,  when  L-  19  mm,  as  a  function  of  portions 
of  a  drop  up  to  an  integral  drop.  The  graph  is  saying  that 
though  in  a  collision,  an  integral  drop  is  expected,  there  is 
still  a  low  level  of  many  fractions  of  drop  fragments  that  may 
be  expected  in  any  firing.  The  high  amplitude  line  for  the 
integral  drop  is  displaced  one  increment  to  the  right  because 
of  plotting  routines. 
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HISTOORAfl  OF  0 .025  FRACTIONS  OF  S  fin  DROPS  INTERCEPTED 

Figure  35.  Water  content  for  fractions  of  5  mm 
Intercepted  drops.  (L  »  19  mm) 

(The  same  general  histogram  occurs  for  other  L  values). 


Figure  36  is  the  associated  cumulative  water  content  for 
Figure  35  (Ls  19  mm).  Again,  normalization  is  used  instead  of 
absolute  units.  All  the  fractions,  from  vei^  minute  to  almost 
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K ISTOORflfl  OF  0 .025  FRACTIONS  OF  S  flfl  DROPS  INTERCEPTED 

Figure  36.  Cumulative  water  content  for  fractions  of 
5  mm  diameter  intercepted  drops  for 
15  mm  diameter  fuze  front.  (L  =  19  mm) 

(Same  general  histogram  occurs  for  other  L  values). 

the  entire  drop,  on  the  average,  contributed  to  almost  half  of 
the  water  content  intercepted  in  any  firing.  About  70',?  of 
water  content  comes  from  drops  intercepted  that  range  from 
3/k  to  a  complete  drop.  Other  values  of  L  spacing  gave  the  same 
general  plot  with  some  small  amplitude  change. 

The  above  approach  is  a  purely  geometrical  analysis  of 
the  encounter.  The  actual  encounter  would  involve  some  smooth¬ 
ing  from  the  interaction  of  the  drops  at  the  periphery  of  the 
fuze  tip  boundary. 


Even  though  in  rsinfields  half  of  the  water  content  (or 
equivalently  half  of  the  impulse  imparted  to  the  fuze)  comes 
from  drops  below  the  median  diameter,  which  is  in  the  vicinity 
of  2  mm,  fuse  response  should  not  be  considered  from  the  aspect 
of  total  momentum  transferred,  but  rather  from  a  subtle  com¬ 
bination  of  the  larger  drops  and  repeated  impacts  that  would 
cause  a  large  amplitude  response  for  those  fuzes  that  are 
more  susceptible  to  a  spring  type  response.  Also  it  is  not 
expected  that  for  each  firing  the  drops  will  be  at  their 
average  spacing  to  produce  a  predictable  resonance  driving 
load. 

A  direct  approach  for  the  facility  would  be  to  design  in 
a  periodicity  or  frequency  important  for  the  type  of  fuse 
considered.  More  important,  whatever  impacts  do  occur,  the 
larger  drops  should  be  used  for  maximum  response  at  any  Impact. 
Five  mm  drops  appear  convenient. 

As  a  specific  case,  consider  once  more  the  example  of  a 
shell  at  3280  ft/sec  in  tropical  rain  at  5.3  in/hr  rate  where 
the  median  diameter  is  2.3  mm.  The  range  is  taken  as  1000  feet, 
and  the  effective  fuse  frontal  diameter  is  15  mm.  The  frequency 
of  impact  is:  Shell  velocity/separation  of  web  stations. 

For  a  shell  at  3280  ft/sec  and  a  50  foot  spacing  of  webs, 
this  is  a  driving  freouency  of  66  HZ,  or  a  time  interval  of 
15  milliseconds.  Engineering  knowledge  of  a  specific  fuze 
would  dictate  what  spacing  to  use  for  the  moveable  web  stations. 
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However,  as  the  expected  frequency  goes  down,  only  a  few 
stations  may  be  possible  along  a  1000  foot  range  to  accomplish 
this  aspect  of  resonance  producing  multi-impacts. 

To  pursue  this  example  further,  the  li3t  below  shows 
how  many  drops  ON  THE  AVERAGE  would  be  intercepted  by  this 
fuze  at  the  cited  drop  diameter  (£D  here  is  0,5  mm). 

Number  intercepted  Frequency  of  impact 
Size  of  drop  ON  AVERAGE  in  1000  (HZ)  for  3280  ft/sec 

(mm)  foot  range  in  and  fuze  diameter  of 

tropic  al  rain.  15  mm. 


1.5 

20 

62 

2.0 

13 

43 

2.5 

7 

23 

3.0 

3.5 

12 

3.5 

1.35 

4.4 

4.0 

0.54 

1.7 

4.5 

0.16 

0.52 

5.0 

0.067 

0.22 

5.5 

0.027 

0.09 

6.0 

0.013 

0.04 

In  this  case  the  driving  frequency  of  the  rain  from  drops 
4mm  and  larger  is  very  low  and  might  not  even  be  pertinent 
for  resonant  or  multi-impact  driving. 

The  drops  in  the  ranre  2  to  3  mm  give  lower  momentum 
transfer  at  each  encounter,  and  if  these  drops  are  uniformly 
spaced,  they  may  be  at  some  important  frequency  of  the  fuze. 
If  all  these  elements  are  present,  then  smaller  drops  may  be 
pertinent.  From  this  uncertainty  comes  the  recommendation 
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for  a  confirmatory  test  at  a  facility  as  Holloman.  This 
issue  may  still  be  resolved  by  a  rain  web  system. 

If  this  frequency  aspect  is  set  aside,  physically  it 
is  a  question  of  whether  many  low  amplitude  impulses,  or  wider 
spaced  higher  ones  are  more  detrimental.  Physically  the 
latter  appears  more  detrimental  from  the  aspect  of  some  minimum 
impulse  needed  to  appreciably  disturb  the  fuae  mechanism. 

Prom  an  applications  aspect,  there  is  the  need  to  drive 
down  the  range  and  replace  the  web  frames  after  each  firing. 

The  lateral  area  of  the  fuze  has  not  been  discussed,  as 
the  momentum  transfer  here  is  degraded  by  the  ogive  contour. 

As  the  ogive  passes  through  the  web,  the  lateral  impacts 
would  be  more  severe  than  in  ordinary  rain  because  of  the 
higher  matrix  density  of  the  web  drops. 


MECHANICAL  SYSTEM  TOR  RAIN  IMPACT  SIMULATION 


For  single  and  even  multiple  rain  impacts,  a  mechanical 
system  such  as  a  rotating  flexible  mechanical  arm  end  weighted, 
or  a  rapid  fire  gun  offer  a  degree  of  simulation. 

There  are  dissimilarities.  Rain  impact  with  high  velocity 
shells  is  an  impulsive  encounter  lasting  no  more  than  about 
ten  microseconds.  Its  erosive  character  is  due  to  high  pressure 
created  at  the  surface  and  by  radial  flow. 

The  surface  pressure  profile  will  not  be  the  same  with  a 
metallic  impact.  There  is  no  radial  flow  and  the  impact 
duration  (depending  on  the  stiffness  of  the  pellet  and  the 
fuze  frontal  configuration)  may  be  an  appreciable  number  of 
microseconds.  This  still  merits  some  consideration.  The  air 
gun  approach  is  straightforward,  but  it  would  be  more  con¬ 
venient  in  a  laboratory  if  the  energy  were  obtained  from  a 
mechanical  rotating  system. 

If  lightweight,  flexible  arms  were  rigidly  attached  to 
a  rotating  shaft,  with  a  spherical  metal  pellet  attached  to 

the  arm  end,  the  following  approach  could  be  used: 

Rotation 


Bending  J 
of  arm  ' 
after 
impact 


u>  -  shaft  angular  velocity 

L  •  distance  from  shaft 
to  fuze 


water  dens  ity 
^ pellet  density 

£  s rebound  coefficient 
1±  6i2 

D  •  rain  drop  diameter 
d  s  pellet  diameter 
?  *  shell  velocity 
▼  *  pellet  velocity  * 
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As  a  first  approximation,  neglect  the  weighs  of  the  flexi¬ 
ble  arm  (this  may  be  a  serious  neglect).  Estimate  the  rotor 
speed,  arm  length  L,  and  pellet  properties  to  simulate  rain 
impact.  First  have  an  equivalence  of  momentum  transfer: 

w«ter  Pellet 

fw!TD3  V  ^»TT  d3 

S  =  5 - 

RPM  a  30**?  (18) 

i r 


Consider  a  specific  case: 


fe  ~  7.6  (steel  pellet)  g/cc 

•  £  s  1.2  (coefficient  of  restitution) 

L  *  0.83  ft  (10  inches) 

D  -  4  nxm  (0.0131  ft)  drop 
V  =  1900  fps 


(RPM)  d3  r  30  ( _1_  \  D3  Vj 
rr  \  7.6  J  &  l 


(RPM)  d3 


30  (0.0131)3  1900 

W  (7.6)  1,2  0.83 


0.00539 


Consider  a 


1000  RPM  shaft: 
d3  *  5.39  x  10~6  (ft3) 


d  -  5.3  a® 


steel  pellet 


Steady  state  operating  conditions  depend  on  the  mounting 
and  stiffness  of  the  rotating  system,  and  experimentation 
is  needed. 
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Duration  of  the  collision  is  obtained  by  monitoring 
electrical  contact  between  the  fuze  (assume  a  metallic  cover) 
and  the  metal  pellet  system. 

If  the  three  pellets  are  so  mounted  as  in  the  sketch, 
the  time  between  collisions  Is: 


20  milliseconds. 


In  the  literature  are  examples  of  water  and  pellet 
type  impact  machines. 
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EXPERIMENTAL- PHOTOGRAPHIC  STUDY  OP  O-C-O  SYSTEM 


Concep t 


i  ;  «  ' 

.  ,  <  '  i 


The  feasibility  of  -in  O-C-O  (overhead-chennel-orifice) 
system  was  briefly  considered,  both  in  theory  and  in  a  simple 

photographic  experiment.  To  produce  a  more 
consistent  drop  distribution,  the  approach 
f ''•  £  ‘yr  of  orifice  flov/  f  or  a  thin  plate  was  used. 

,  .  (  i 1  1  Water  is  in  an  overhead  channel  at  a 

!  i  '  i  / 

ill  1  7  liquid  level  with  holes  of  various  sizes 

j  am  i  I  i 

STiy  drilled  into  the  channel. 

4  '  Questions: 

'What  si  ze  holes  should  be  used? 

How  uniform  are  the  emerging  drops? 

What  height  is  needed  ford  rop  formation? 

What  is  the  drop  spacial  separation? 

Is  such  a  system  maintainable? 

The  veloci ty  at  vh  1  ch  water  emerges  as  a  stream  from  a 

/  /  //N)/  ■  /  /  |  b  fluid  height,  h,  is  close  to  ideal  gravity  fall 

*  v  v  it  0.99  V  2ghn 


The  amount  (discharge  rate  Q)  of  water  depends  on  the 
orifice  configuration,  as  the  effective  exit  area  is  not  the 
geometrical  area.  Por  a  thin  plate  and  straight  orifice,  the 
discharge  coefficient,  f,  is  cited  about  0.6  (dimensionless). 

The  discharge  rate  Q  (ft^/sec)  through  the  orifice  is: 

Qsf  AV~2gh~1  ,  where  A  =.  geometric  orifice  area. 

Consider  that  at  some  distance  down  from  the  orifice,  a 
stab]e  drop  pattern  has  been  produced.  Assume  all  the  dm  ps 
are  of  the  same  size,  and  are  equally  spaced.  These  are  poor 
assumptions,  but  give  an  initial  insight. 


The  drop  diameter  D  Is  not  expected  to  be  the  same  as 


the  orifice  diameter  d.  Some  function  should  exist,  and 


consider  thi3  as  linear  with  C  as  the  slope: 

D=  C  d.  At  some  region  below  the  channel,  a 
steady  state  exists.  Mass  rate  balance  between  the 
number  of  drops,  n,  passing  this  location  per  unit 
time,  and  the  efflux  from  the  orifice  is  the  same: 
n  TT d3/6  =  f  A  =  f 'TT d2  /U  (19) 

The  number  of  drops/sec  passing  this  region  becomes: 


n  -  3  f V2  g  h 

2  C2  D 

Or  conversely,  the  discharge  coefficient  is:  f-  2  n  ^  ^ 

3  (2TT* 

Single  flash  photography: 

'  #  This  could  obtain  the  coefficient  C  and  n: 

•  L°  length  of  section  in  stabilized  drop  formation 
N=  number  of  drops  in  this  distance  L 
^  «  v-=  velocity  of  drops  in  this  section  L 

Number  of  drops  passing  per  unit  time: 
nsN/At  =  N  V/L 

^  v  Inspection  of  the  drop  size  should  give  the  slope  C. 
Of  concern  is  the  spacing  between  the  drops,  which  will 
be  too  close :  3  = separation  of  the  drops: 


s  -  h  -  L  L 


2  v  C  D 

3  f^gh1 


For:  v=l5  fns  )  33I4..  23  D 

<3*1.2  - 

f-  0.6 

h  -  6  Inches 


D(mm)  S(centlmeters  \ 

2  0.85 

4  1.7 


The  conclusion  is  that  unless  the  drops  are  scattered 
horizontally,  thoy  are  too  closely  spaced  vertically.  Their 
speed  is  too  slow. 


Experiment  and  Tabulation  of  Results 


A  sample  photographic  setup  to  observe  drop  formation 
from  a  channel-orifice  system  consisted  of  four  large  cans, 
with  carefully  drilled  orifices.  Five  feet  below  the  bottom 
of  the  can,  a  35mm  flash  camera  observed  the  drop  formation 
against  a  black  background. 

Skim  milk  had  to  be  used  for  clear  identification  of  the 
drops.  The  flash  duration  was  less  than  20  microseconds 
(no  visible  distortion  from  motion  could  be  seen).  Appendix  A 
has  several  samples  of  the  60  SINGLE  flash  pictures  taken. 


5  feet 


Gan 

1 

2 

3 

4 


Four  different  cans  were  used: 

Number  Hole  diameter 

of  holes  inch  inch  mm 


2 

2 

2 

1 


1/16 

3/32 

i/e 

3/16 


.063 

.094 

.125 

.188 


1.6 

2.38 

3.1 

4.76 


The  liquid  level,  h,  was  maintained 
at  either  3  of  6  inches,  representing 
exit  velocities  of  4  and  5.7  fps .  The 


experiment  was  dene  under  quiet  condi¬ 
tions  in  the  fluid.  With  rotational 
motion  in  the  fluid,  the  stream  was 
scattered  sideways  in  what  appeared  a  random  manner.  Unfor¬ 
tunately,  no  pictures  ’were  taken  of  this  situation. 
Experimental  results: 

1.  An  upper  limit  liquid  drop  size  is  created  by  the 
orifice.  This  appeared  independent  of  the  liquid  level  used. 

2,  A  sharp  peaked  blmodal  drop  size  distribution  was 
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9 

0 

0 

0 

o 


created.  This  could  be  seen  in  a  picture  of  the  stream 
about  a  foot  belovv'  the  orifice  where  the  drop  for¬ 
mation  was  essentially  completed.  The  larger  sine 
drop  is  larger  than  the  orifice  diameter. 

3.  Drops  are  not  in  uniform  vertical  formation, 
and  are  too  closely  spaced  vertically  in  some  pic  tures. 

Some  swirling  motion  in  the  fluid  might  be  desireable 
for  horizontal  spreading  of  the  drops. 

4.  Drops  were  photographed  in  a  vertical  distance 

of  7  inches,  at  a  location  5  feet  from  the  orifices.  The  drops 
from  the  largest  orifice  were  not  very  spherical,  and  it  was 
difficult  to  estimate  their  diameter. 

'  Table  12  lists  the  holes,  drop  diameters  (measured  with 
a  ruler  from  8x10"  pictures),  and  the  average  number  and 
standard  deviation  of  the  number  of  drops  seen  over  these 
7  inch  intervals.  These  are  averages  for  all  experiments, 
both  6  and  3  inch  channel  fluid  height. 


Tahle  12.  Drop  sizes,  number,  liquid  content  from 
orifice  flow  experiment 


Bole 
(mm  ) 

Dron  diamet< 

1  Large  Med 

I 

ar  (mm) 

Small . 

Numb  9i 
7"  drc 
below 
Large 

1  .drc ps  < 
>P,  5  fei 
orifice! 
Med 

..  —  ,  -j 

aver 

9t 

3 

Small 

'Orai 

(G) 

3"-6" 

ms 

6" 

i 

SO 

>  SO 

-5  so 

1  1.6 

2  2.38 

3  3.1 

4  4.76 

2.27  :  — 
3.18  '  — 
3.52  ■'  — 

*5.5  -3.6 

mm 

n 

9.1  2.0 

7.1  1.4 

7.2  1.3 

6.5  2.0 

1 

3.3  1.4 

8.3  1.6 

5.4  1.3 

5.5  1.0 
3.8  1.5 

.058 

.129 

.182 

-.656 

.066 

.135 

.184 

-.713 
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The  last  two  columns  In  table  12  give  the  average 
number  of  GR4M3  of  liquid  seen  in  this  7  inch  interval  for 
two  conditions: 

(1)  For  all  the  data,  both  3  and  6  inch  heights.  Forthe 
four  type  orifices,  the  meai  fluid  height  was: 

Orifice  Mean  fluid  height  (in) 

1  TTF 

2  5.3 

3  5.3 

k  3.8 

(2)  Data  for  only  the  6  inch  fluid  height. 

To  test  the  consistency  of  the  data,  the  discharge 
coefficient,  f,  was  calculated  for  the  orifices,  which 
theoretically  should  be  0.6.  Since  the  drop  velocity  could 
not  be  determined  experimentally,  a  velocity  of  17  fps  was 
used  for  these  drops.  iee  the  section  on  Reynolds  number 
and  figure  3. 

(1)  Calculate  liquid  flow  through  seven  inch  window: 


R  (cc 


7777)1  =f  G  ( grr  m3  )1  f  17  ftlfia  inT  [— ~ 

1  17"  wine  ow  JL  sec  JL  ft  J 


(grams  )\  (21) 


(2)  Calculate  the  liquid  flow  from  the  orifices: 

Q  (cc/sec)  -  f  A  ^2  gif  -  ^  1|  fcif.v  *1  JL/jr'  li-- ^  ^ 

Q  (cc/sec)  -  1240  d2(ln)  fh  (ftf 


d  -  orlfi  ce  diameter  l  in) 
b -fluid  heigiit  in 
channel  (ft) 

f -discharge  coefficient 


Find  f: 

f  -  29  0 
W 


The  coefficient  !.V  is  calculated  for  the  various  conditions. 
Mean  height  IT  is  used  for  the  combined  3  and  6  inch  heights. 
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Hole 

d(in) 

W 

(&> 

W 

(h  =  0.5) 

CD 

(b) 

0 

(h  =  0.5  ft) 

1 

O.O63 

3.11 

3.48 

0.54 

0.55 

2 

0.094 

7.28 

7.75 

0.51 

0.51 

3 

0.125 

12.9 

13.7 

0.41 

0.39 

4 

0.188 

24.7 

31. 

0.77 

0.67 

There  is  crude  agreement  with  the  theoretical  value 
f=0.6.  The  value  is  lower  for  the  third  hole;  and  it  was 
difficult  to  estimate  the  drop  diameter  for  the  fourth 
size  hole. 

Application  to  outdoor  test 

1.  A  spec’flc  drop  size  distribution  can  be  created. 
The  rain  rate  can  be  controlled  to  a  degree,  and  is  propor¬ 
tional  to  the  square  root  of  fluid  height  in  the  channel. 

2.  The  vertical  drop  separation  is  too  close  (at  times 
the  drops  are  almost  binary  and  would  be  equivalent  to  a 
very  large  drop  on  a  fuze).  Some  swirling  motion  in  the 
channel  may  change  this. 

3.  To  prevent  clogging  of  the  orifices,  a  covered 
and  filtered  channel  system  v.'ould  be  needed. 

4.  For  an  inten 3e ,  calibrated  rainfield  of  short 
length,  this  O-C-O  system  may  be  a:  i table. 
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APPENDIX  A.  Photos  from  study  of  0-0-0  system. 

Of  the  35  xty  photos  taken  for  the  overheed-channel- 
orifice  system,  only  five  are  Included  here.  The  liquid 
height  in  the  channel  for  these  photos  was  six  inches. 

A  centimeter  ruler  is  referenced  with  each  photo. 


Figure 

Number 
on  photo 

Hole  size 
number 

Hole  Large  drop 

Orifice  (  am)  d lame ter (mm) 

37 

1 

1 

1.6  (note  "binary" drops )  2.27 

38 

3o 

1 

1*6  2.27 

39 

59 

2 

2.38  3.18 

bo 

22 

3 

3*1  3.52 

41 

39 

b 

4*76  (drops  unstable)  '£.(5.5) 

8  0 


Figure  37.  Drops  from  1.6  mm  orifice 


(Large  drop  diameter  approximately  2.27  ran) 
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Figure  38*  Drops  from  1.6  mm  orifice. 

(Large  drop  diameter  approximately  2.27  ram) 


82 


Figure  39. 


Drops  from  2.38  ram  orifice. 

(Large  drrp  diameter  approximately 
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figure  liO. 


Drops  from  3.10  mm  orifice. 

(Large  drop  diameter  approximately  3»52  mm) 
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Figure  ip..  Drops  from  4*76  mm  orifice, 

(Large  drop  diameter  approximately  5.5  ms) 
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RAINDROP  IMPACT  ON  PD  FUZES 
(UNCLASSIFIED) 


by 


J.  K.  Domen 

Picatinny  Arsenal 
Dover,  New  Jersey  0780 1 


ABSTRACT 


A  summary  is  given  of  analytical  expressions  related  to 
rainfields,  and  of  specific  information  (graphs,  tables)  of  reported 
rain  drop  size  distributions  per  unit  volume  for  various  natural  rains 
and  artificial  rains  at  testing  facilities.  A  review  is  made  of  rather 
diverging  analytical  approaches  to  ascertain  the  pressure  and  force  time 
pulse  imparted  by  drops  on  rather  riqid  surfaces. 

A  computer  simulation  which  employed  momentum  conservation  is 
presented  for  a  specific  point  detonating  artillery  fuze  nose  (spring- 
mass)  response  to  hiqh  velocity  flight  in  heavy  natural  and  excessive 
artificial  rainfall,  with  corresponding  experimental  results  at  Hollo¬ 
man  Air  Force  Base. 

A  pertinent  summary  is  given  of  experimental  results  of  erosion 
of  rain  impacted  surfaces,  the  basic  phenomenon  involved,  methods  employed 
for  erosion  abatement,  and  analytical  considerations  for  erosion  process. 


RAINFIEID  CHARACTERISTICS 


For  both  naturally  occurring  and  artificially  produced  rains,  raindrop  diameters 
range  from  a  fraction  of  a  millimeter  to  rarely  over  six  millimeters.  The  liquid  water 
content  L  (g/m3)  of  rain  is  defined  as  water  mass  M  for  some  given  volume  V  where  the  sum 
is  taken  over  all  drops  in  the  volume: 

l  -  M/v  =  HE—  £  Nj  Dj3  (g/m3)  (D 

6V  <•/ 

where  F  -  Final  group  of  drops. 

Nj  •  Number  of  drops  with  average  diameter  Dj  per  volume  V. 

P  •  Water  mass  density  (I  g/cc). 
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The  rainfall  accumulation  rate  R  (in/hr)  is  then  given  as: 


t&n.  i  a: 


R  -  x/t  =  — IL —  r  w;  N;  (mm/hr) 

6V 


where  t  *  Time 

v;  »  Terminal  velocity  of  drop  of  diameter  0 
x  -  Accumulation  depth 


(2) 


For  quiet  conditions,  all  drops  attain  their  terminal  velocity^  after  a  fall  of 

about  50  feet,  and  this  speed  can  be  expressed  approximately  as: 

3. - - 

v  (ft/sec)  *  16  V  D(mm)  +  2  (3) 

for  drops  in  the  diameter  range  of  I  to  6  millimeters.  A  droplet  of  1.5  mm  diameter  requires 
about  20  foot  free  fall  to  attain  an  18  ft/sec  terminal  velocity.  Table  1  lists  the  reported 
terminal  velocity  of  drops  . 

TABLE  I.  WATER  DROP  TERHINAL  VELOCITY 


0(mm) 

v  (ft/sec) 

D(mm) 

v  (ft/sec) 

1.25 

15-9 

3-25 

27-2 

1.50 

18.1 

3-50 

28. 

1.75 

20. 

3-75 

28.6 

2.00 

21  .6 

4.00 

29. 

2.25 

23. 

4.50 

29.8 

2.50 

24. 3 

5.00 

30.3 

2.75 

25.5 

5-55 

30.5 

3-00 

26.4 

6.00 

30.5 

Water  drops  in  an  artificial  field  generally  have  a  lower  terminal  velocity 
because  of  the  small  drop  height,  and  consequently  a  higher  drop  concentration.  An 
equ i valent  rate  for  the  artificial  field  is  defined  as  the  rain  rate  obtained  if  the  drops 
were  travelling  at  their  terminal  velocity,  ana  is  higher  than  the  usual  accumulation  rate. 

The  mean  diameter  ranges  from  1.4  to  2 . I  nm  for  natural  rain  and  is  that  diameter 
drop  having  a  mass  equal  to  the  mass  content  M  divided  by  the  total  number  of  drops  N  per 
unit  volume:. 


where  N  *  LNj  «  Total  number  of  drops  in  volume  V. 


68 


The  med  i  an  d  i  ano  ter  ranges  from  ]  .2  to  2  *nr  for  natural  rain  and  is  that  diameter 


0.864 

l(g/m3)  =  0.075  R  (mm/hr 


Logarithmic  least  squares  fit  for  Miami,  Florida  data  has  been  reported  as: 


(5) 


drop  for  which  half  the  total  mass  M  :s  above  this  diameter  and  half  below.  It  occurs  when 
during  the  Summation  from  either  spectrum  end  the  equality  occurs: 

f  t 

L  N.D.3  =  %  t  NjD.3  (D— -M...) 

i=1  \ 

An  empirical  logarithmic  relation  appears  to  exist  between  rainfall  rate  and 
liquid  content . ^  Data  for  the  range  from  moderate  rain  (R  3  mm/hr)  to  cloudburst 
(R  *?  100  mm/hr)  can  be  expressed  as. 


(6) 


.(4) 


0.95 

L(g/m3)  =  0.052b  p  (mm/hr)  (7) 

Natural  rain  data  indicates  that  tne  median  or  maximum  drop  diameter  (mm)  is  pro- 

(M 

portional  to  the  cube  root  of  tne  rain  rate  (mr./i  r) 

Natural  rain  spectrum  of  frequency  of  occurrence  of  drop  diameter  usually  has  a 
peak  between  0.7  and  2  mm.  If  the  frequency  of  drop  size  for  natural  rain  is  plotted  versus 
tne  log  of  the  drop  diameter,  the  curve  appears  Gaussian  normal  (Figure  1).^  Natural  rain- 
fall  appears  to  plot  well  with  the  express ioi 


Nq  dD  = 


N 


e 


In  (D/G) 
In  S 


2 


d  (In  D) 


(8) 


where  M  *  Total  number  drops  per  unit  volume  (area  under  either  curve)  ■  ;.Nj 
Mr,  r  Number  drops  per  unit  volume  per  dirnneter  interval 
N|  “  Number  drops  in  diameter  irter/.il  hi  centered  at  diameter  Dj 


rt9 


V«T  *?»■*•■ 


G  -  Geometric  mean  diameter 

in  G 


S  »  Geometric  standard  deviation 

(In  S)2  -  I.  C  Nj  (In  0,)  -(in  G)2 


G  -  N[  In 

N  4i 
» 

ation 

-  i  ^Nj  (In 


For  a  cited  distribution  of  114  flm/hr  at  the  Canal  Zone,  these  values  are: 


In  G  -  .3414 

G  -  I .407  mm 
Mean  Diameter  ■  1.81  mm 


In  S  -  .414 

S  -  I.5U 

Median  Diameter  -  2.4  mm 


txcept  For  small  diameter  drops,  another  Form  has  been  used  For  Nq  (drops  per  cubic 
centimeter  per  centimeter  diameter)  but  has  had  limited  success:^ 

N„  -  0.08  cm-4  e  ■*° 


where  X(cm_1)  -  41  R  (mti/hr)  . 


Natural  RainFall 


Rain  drops  usually  occur  below  20,000  Feet  but  may  be  Found  as  hiqh  as  50,000  Feet. 
They  seldom  appear  above  35,000  Feet.^)  Drop  size  distributions  are  ordinarily  reported  and 
plotted  In  histoqram  Form  as  the  number  oF  drops  Nj  per  volume  V  within  a  certain  diameter 
interval  AD  near  a  diameter  Dj  versus  diameter  D.  A  "standard"  rainFall  (Table  2)  has  been 
reported  with  90$  conFidence  that  In  95$  of  the  area-month  combinations  of  the  U.S.,  99$  of 
all  rainFall  intensities  will  be  equal  to  or  less  than  this  standard  oF  2  in/hr. ^ 

TABLE  2.  DROP  SIZE  DISTRIBUTION  OF  "STANDARD"  RAIN  (AO  -  .25  mm) 


D| (mm) 

Nj  (Drops/m3) 

D|  (mm) 

N  |  ( Droi 

.125 

0 

3.375 

5.13 

.375 

1026 

3.625 

3.07 

.625 

452 

3.875 

1 .84 

.875 

227 

4.125 

1.23 

1.125 

145 

4.375 

0.67 

1.375 

107 

4.625 

0.41 

1.625 

76 

4.875 

0.25 

1.875 

5° 

5.125 

0.17 

2.125 

43 

5.375 

0.11 

TABLE  2.  DROP  SIZE  DISTRIBUTION  OF  "STANOARD"  RAIN  (AD  -  .25  tin)  -  CONT'D 


D| (mm) 

Nj  (Drops/m3) 

Dj  (mm) 

N  |  (Drops/m3) 

2.375 

29 

5.625 

0.08 

2.625 

20 

5.875 

0.04 

2.875 

12 

6.125 

0.02 

3.125 

7-69 

6.375 

0.02 

N  -  2207.7 

In  Fiqure  2  is  summarized  for  various  world  locations,  the  rainfall  rate  equalled  or 
exceeded  for  an  indicated  percentaqe  of  total  time  In  that  location,  not  just  during  the  time 

(q) 

it  is  actually  raininq.  ■“  This  data  is  based  on  approximately  one  year,  exceot  Panama  Is  ex- 
aqqerated  as  data  here  is  taken  durinq  the  wet  season  from  June  to  November.  Table  3^  lists 
the  rates  not  for  Just  the  percentaqe  of  total  time  (Fiqure  2)  but  also  for  the  percent  of  time 
given  that  It  Is  raininq.  Table  3  includes  rates  for  desert  and  temperate-tropical  regions. 
Specific  data  of  natural  rain  drop  sizes  for  the  Canal  Zone  (Table  4  and  Fiqure  3)^'®',  Miami 
Table  5)^''^,  Marshall  Islands  (Table  6) ^  and  North  Carolina  (Fiqure  4)^  was  obtained  by 
continuously  photograph i nq  a  small  volume  and  obtalnino  the  number  and  drop  size  per  volume 
(m3) .  The  aver age  distributions  for  all  major  climatic  zones  are  made  from  data  at  nine  loca¬ 
tions  around  the  world  with  rates  within  +_  127,  of  the  average  rate,  and  are  plotted  in  Fiq¬ 
ure  5.  Here,  927  samples  were  used  for  the  5.2  rnn/hr  rate  and  154  for  the  95.6  im/hr. 

TABLE  3.  RAINFALL  RATE  (rmt/hr)  EQUALLED  OR  EXCEEDED  FOR  INDICATED  PERCENTARE  OF 


Spaclal  Distribution  of  Raindrops 


Quantitative  information  is  not  readily  available  for  the  spacial  distribution 
between  drops.  Photographs  do  Indicate  non-uniform  spacing  and  for  simulation  purposes  a 
Poisson  distribution  might  be  reasonably  assumed.  For  example,  for  a  PD  Fuze  with  a  nose  of 
effective  frontal  normal  area  A  at  high  velocity  v  through  a  rainfield  of  N  drops  per  volume 
V,  consideration  of  the  cylindrical  volume  swept  out  by  the  nose  gives  the  average  time 
between  collisions  (TBCOL) : 


T8C0L  -  1 / (NAv)  (10) 

and  the  average  distance  between  drops  for  a  uniform  drop  spacing  is  1 / { NA) .  If  N  *  10,000 
drops/m3,  and  the  nose  is  0.5"  in  diameter,  the  average  distance  becomes  0.787  m  -  2.58  ft/per 
impact.  With  the  Poisson  term  NA  «  0.387  fraction  drop  impacted/foot  travel,  the  probability 
of  M  impacts  in  one  foot  of  travel  becomes: 

e  ~NA  (NA)  M 


P(M)= 


(ID 


M! 


P(0)  -  .6791 
P( I)  -  .2628 

P(2)  *  .05085  (Probability  of  2  drops  in  one  foot  travel  «  5.U) 
P ( 3)  -  .00656 


Artificial  Rainfall 

Data  is  reported  for  two  facilities  at  Itol  loman  Air  Force  Base,  New  Mexico: 

(I)  An  artillery  range  where  H  -  1/2  U  80200  (deluge)  nozzles  are  usually  employed  at  3  1/2 
ps  i  with  each  nozzle  discharging  5.9  gal/min.  The  pipe  supported  nozzles  stand  on  the  ground, 
alternating  on  opposite  sides,  and  spray  at  about  65  degrees  above  the  horizontal.  The 
accumulation  rate  from  the  nozzles  is  from  about  20  to  36  in/hr.  The  drop  intensity  Is 
nine  times  as  heavy  as  tropical  rain  of  1 32  mm/hr.  (2)  Rocket  mono-rail  sled  track.  35,000 
feet  long  of  which  6000  feet  is  through  a  rainfield  produced  ordinarily  by  the  standard 
nozzle  H  -  l/A  U  8070  which  gives  a  smaller  number  of  drops  and  fewer  large  drops.  Sled 
velocity  up  to  Mach  5  is  possible.  If  the  deluge  nozzles  are  employed,  only  about  2000  feet 
of  the  track  can  be  used  because  of  the  high  delivery  rate. 
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The  nozzles  produce  an  excessive  number  of  small  drops  (about  0.5  tnn  diameter)  and 
there  Is  Insufficient  drop  heiqht  (especially  for  smaller  size)  to  attain  terminal  velocity. 

The  average  drop  size  spectra  for  these  facilities  with  their  standard  nozzles  are  listed  in 
Table  7-^  Other  reported  "deluge"  nozzle  (H  1/2  U80200)  distributions  are  in  Table 
The  average  standard  artillery  ranqe  and  rocket  track  rainflelds  are  plotted  in  Fiqure  6  and 
compared  with  tropical  rain  of  132  mm/hr.  The  natural  rain  falls  between  the  two  distributions 
for  diameters  greater  than  2.6  mm.  Figure  shows  scatter  for  the  rocket  track.  More  Informa¬ 
tion  is  contained  in  reference  9. 

(tit) 

Table  9  Is  a  Sandia  distribution.  ' 

FIGURE  2.  THE  FREQUENCIES  OF  OCCURRENCES  OF  ONE-MINUTE  RAINFALL  RATES 
AT  VARIOUS  LOCATIONS  III  TERMS  OF  PERCENTAGE  OF  TOTAL  TIME 
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FIGURE  3.  AVERAGE  RAINDROP  DISTRIBUTIONS  FOR 
THE  PINA  RANGE,  CANAL  ZONE,  JUNE  2 


3 
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FIGURE  It.  AVERAGE  DROP  SIZE  DISTRIBUTIONS  FROM  NORTH  CAROLINA 
FOR  VARIOUS  RAINFALL  RATES,  HM/HR 


FIGURE  5.  AVERAGE  DROP-SIZE  SPECTRA  FOR  NATURAL  RAINFALL  RATES 
OCCURRING  O.OI,  0.1,  0.5  and  1.0*  OF  THE  TIME 
(These  curves  are  for  all  available  data  from  nine 
locations  In  all  major  climatic  zones  around  the 
world.) 
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NUMBER  OF  DROPS/W5  C.l 


FIGURE 


i 


7.  THE  AVERAGE  SPECTRUM  FOR  THE  STANDARD  ROCKET  TEST 
TRACK  RAINFIELD  BASED  ON  FOUR  NO-WIND  TEST  RUNS 
(The  vertical  bars  show  the  ranqe  of  measurements 
from  the  four  runs.) 


i 

1 
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TABLE  4.(C0NT'D)  AVERAGE  RAINDROP  DISTRIBUTIONS 
CANAL  ZONE  (JUNE -JULY  1968) 

(AO  "  0.2mm) 


R(mm/hr)  Rate 
L(g/m3)Liquid  Content 
DH(nn)  Median  Diameter 


NS  Number  m  3  Samples 
N  Drops /m 3 


R 

11.7 

18.8 

30.5 

49.4 

77-9 

114. 

L 

■  57 

.86 

1.35 

2.1 

33 

4.86 

DM 

1-9 

2.1 

2.2 

2.4 

2.5 

2.4 

NS 

145 . 

82. 

73. 

64. 

59- 

15. 

N 

22<t.5 

274.9 

380.1 

538.4 

883-9 

1461 .2 

Dj  (mm) 

•  55 

•  76 

1.12 

2. 

4.9 

14.2 

24.7 

■  75 

A.  3 

7.15 

11  . 

25.9 

60.2 

124.2 

■  95 

18.2 

20.3 

30.2 

55-5 

120.6 

247.4 

1.15 

46.8 

69.2 

124.1 

226.3 

1.35 

57.75 

69-3 

104.6 

181  . 

1.55 

mm 

fill 

58.2 

65.2 

85.6 

136.2 

1.75 

31-9 

40.8 

51.6 

62.05 

76.9 

1.95 

19.7 

29.4 

42.9 

55-1 

77-3 

2.15 

11.7 

19-5 

28.7 

43-06 

65.9 

90. 

2.35 

6.03 

11.3 

18.6 

29-3 

68.6 

2.55 

2.86 

6.8 

12.6 

19-4 

ii^plSfeiin 

49.2 

2.75 

1.62 

4.1 

7-0*4 

12.6 

u 

32-4 

2.95 

.88 

2.43 

5. 

8.9 

15.7 

24.3 

3.15 

.52 

1.17 

3.17 

6. 

10.5 

14.5 

3- 35 

•32 

.68 

1.5 

4.16 

7.1 

10.2 

3-55 

.11 

.45 

1.35 

3-14 

4.5 

7.6 

3-  75 

.08 

.24 

■  7 

1 .67 

2.71 

4. 

3-95 

.04 

.22 

.44 

1 .04 

1.35 

2.68 

A. 15 

.02 

•  09 

.24 

•  91 

.89 

1  .4 

A. 35 

.02 

.04 

.13 

•  34 

.62 

.84 

A. 55 

.02 

.04 

.24 

.27 

.49 

A.  75 

.01 

.04 

.2 

.25 

.21 

A. 95 

.04 

•  13 

.09 

.14 

5.15 

.06 

.2 

.07 

5-35 

.01 

•  05 

.14 

5- 55 

.01 

.02 

.2 

.07 

5.75 

.01 

.02 

£.35 

.01 

.01 

.02 

.07 

6.75 

.02 

100 


TH  «  Thunderstorm 


TABLE  5.  RAINDROP  SI  It  DISTRIBUTION 
Ml  AMI ,  FLORIDA 

CR  •  Continuous  Rain  RS  “  Rain  Showers 


R  *  Rain 


Wet  Season 


Dry  Season 


D(mm) 

TH 

CR 

RS 

TH 

RS 

R 

•  5 

185 

0 

2 

85 

0 

0 

.6 

220 

4 

0 

252 

1 

0 

•  7 

387 

12 

0 

673 

4 

0 

.8 

589 

29 

1 

1032 

6 

2 

■  9 

1029 

101 

6 

1883 

6 

4 

1.0 

1393 

210 

41 

1426 

13 

29 

1.1 

1670 

365 

124 

1843 

25 

102 

1.2 

1966 

407 

238 

1778 

40 

171 

1.3 

2087 

516 

561 

2121 

41 

292 

1.4 

21 1 1 

612 

1107 

2188 

59 

407 

1.5 

2072 

736 

1837 

2303 

68 

454 

1.6 

2081 

803 

2482 

2223 

102 

429 

i  •  7 

611 

2316 

1.8 

1.9 

1517 

19  1. 

2.0 

1280 

:  MB 

2.1 

999 

309 

719 

1411 

79 

348 

2.2 

891 

202 

460 

933 

37 

225 

2-3 

696 

191 

399 

767 

41 

194 

2. A 

658 

159 

325 

586 

25 

154 

2.5 

540 

157 

237 

557 

28 

100 

2.6 

433 

159 

192 

516 

16 

63 

2.7 

342 

108 

104 

441 

23 

55 

2.8 

296 

113 

52 

215 

18 

47 

2.9 

281 

65 

41 

172 

20 

25 

3-0 

218 

64 

24 

107 

II 

16 

3.1 

187 

47 

7 

too 

6 

13 

3-2 

163 

32 

2 

87 

8 

14 

3.3 

101 

14 

3 

66 

5 

7 

3.4 

98 

17 

1 

41 

2 

5 

3.5 

87 

19 

2 

50 

2 

4 

3-6 

96 

9 

28 

1 

1 

3.7 

55 

6 

25 

0 

4 

3-8 

68 

7 

12 

2 

4 

3.9 

64 

3 

8 

0 

4 

4.0 

49 

1 

16 

2 

2 

4.1 

50 

4 

17 

1 

4.2 

39 

1 

7 

1 

4-3 

23 

1 

1 1 

4.4 

19 

2 

5 

4.5 

13 

0 

6 

4.6 

15 

0 

4 

4.7 

5 

0 

3 

4.8 

10 

1 

4 

4.9 

7 

1 

5.0 

9 

2 

5.1 

3 

2 

5.2 

7 

2 

5.3 

6 

0 

5.5 

4 

2 

5.6 

3 

2 

5.7 

1 
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TABLE  6.  AVERAGE  RAINDROP  SIZE  DISTRIBUTION 
MARSHALL  ISLANDS 


N(/m3) 


)  1 

2 

6 

7-1 

12.6 

26.8 

67-3 

96.8 

170.5 

9^.87 

169.56 

252.6 

370.33 

615.32 

869.66 

1212.66 

2120.22 

2620.17 

5.61 

8.61 

6.03 

2.59 

7.77 

6.6 

16.67 

128.53 

25. 

7.29 

12.37 

8.07 

7.06 

21  .32 

13.12 

37.11 

90.3 

63.13 

10.16 

17-33 

15.35 

16.17 

33.62 

21  .96 

56.96 

167-59 

59.62 

12.76 

21.33 

26.56 

32.08 

69.69 

61  .18 

79.73 

158.6 

120.68 

13-32 

22. 

28.66 

36.95 

69.6 

62.13 

81.88 

195.89 

125-86 

15.02 

26.25 

65.98 

62.17 

85. 

99.85 

113.62 

I83.8I 

233.6 

10.69 

19.36 

35.59 

56-51 

79.52 

106.27 

96.66 

135.65 

222.65 

7.22 

I6.63 

28. 

66.52 

70.95 

98.3 

95.26 

III .69 

187.22 

6.67 

9.22 

19.71 

36.78 

58.5 

97-18 

91.66 

106.52 

152.5 

2.88 

6.17 

13.23 

26.23 

68.77 

83.89 

87.91 

96.95 

133-8 

1.92 

3-97 

9.77 

16.85 

36.65 

68.62 

79.1 

98.15 

112.52 

1.15 

2.55 

6.65 

11.3 

25-2 

53.61 

76.56 

90.3 

119.85 

•  76 

1 .71 

6.05 

7.6 

17.33 

62.1 

63.91 

86.17 

118. 71 

.56 

1.37 

2.95 

5.69 

11.37 

35.1 

59-38 

93-3 

123.67 

•  31 

.88 

1.87 

3-52 

7.75 

20.38 

65.55 

65.09 

82.65 

.26 

.65 

1.29 

2.62 

5.51 

16.06 

37-8 

72.12 

96.91 

.13 

•  39 

.78 

1 .82 

3-12 

8.81 

25-83 

65-25 

72.53 

•  07 

.29 

.57 

1.08 

1.95 

6.29 

21  .32 

62.15 

69.08 

■  05 

.22 

.63 

.83 

1.25 

3.66 

13.69 

30.79 

53-73 

.05 

.16 

.6 

.63 

•  78 

2.55 

10.16 

30.79 

56.72 

.02 

.07 

.23 

.51 

.62 

1  .61 

6.63 

20.77 

38.33 

.02 

.07 

.12 

•  38 

■  67 

1. 

5.11 

21.39 

37-92 

.01 

.06 

.09 

.19 

.28 

■  37 

3-19 

16.88 

29.86 

.01 

.03 

.05 

.23 

•  29 

.52 

2.81 

9-71 

22.96 

.01 

.02 

.03 

.11 

.25 

.26 

1 .66 

7.66 

18.29 

.01 

.05 

.05 

.18 

•  03 

1 .21 

6.3 

15.7 

3 

.01 

.07 

0 

.03 

.07 

1 

.01 

.01 

1 

.01 

.00 

1 

.01 

.02 

.00 

.01 

6 

1.18 

3.2 

10.02 

9 

.86 

2.58 

8.67 

3 

.57 

2.17 

6.82 

9 

.23 

1.03 

3.62 

5.3  1 

.1 

.1 

5-6 

.21 

.1 

5-5 

.01 

.00 

.00 

5-6 

.00 

.00 

.1 

5-7 

.01 

.21 

TABLE  7-  HOLLOMAN  AVERAGE  DROP  SIZE  SPECTRUM 
ARTILLERY  RANGE  (STANDARD  NOZZLE  H-l/2  U  80200) 


D(inm) 

No./m3 

0.1  m 

No.  Drops/m3 
<_  D 

Water  In  Drops  < 
(g/m3) 

D 

0.5 

6059 

6059 

•  36 

0.7 

1650 

7586 

•  98 

1.0 

820 

11089 

2.51 

1.5 

360 

13663 

5.69 

2.0 

210 

16763 

9.70 

2.5 

160 

15573 

15.28 

3-0 

_ 

16150 

22.18 

3.5 

60 

16517 

29.32 

6.0 

33 

16726 

35.67 

^-5 

17 

16860 

60.29 

5.0 

7-9 

16897 

63.68 

5.5 

2.6 

16922 

65.60 

6.0 

1.6 

16932 

66.65 

6.5 

0.5 

16936 

67.12 

All  Sizes 

16961 

67.68 

ROCKET  TEST  TRACK 

(STANDARD  NOZZLE  H-l/6 

U  8070) 

(NEGLECT  DROP 

<  0.5MM  IN  DIAMETER) 

No.  Drops/m3 

Water  in  Drops  < 

D| 

D(mm) 

<_  D 

(g/m3) 

0.75 

1650 

0.209 

1.25 

3110 

1  .256 

1-75 

6365 

3-617 

2.25 

6882 

5-938 

2.75 

5032 

7. '65 

3-25 

5068.9 

7-687 

3-75 

5076.6 

7.860 

6.25 

5078.6 

7-927 

6.75 

5079.2 

7.958 

Artillery  Range 

Rocket  Track 

L(g/m3) 

67-7 

7.96 

Equivalent  Rainfall  Rate 

1257mm  (69.5  in)/hr 

158.7mm  (6.25  !n)/hr 

Mean  Volume  Diameter  (mm) 

1  .75 

1  .66 

Median  Volume  Diameter  (nri) 

3-21 

1.8 
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TABLE  8.  HOLLOMAN  DROP  SIZE  SPECTRUM 


(H  1/2  U  80200  NOZZLES  AT  6  PSI  MANIFOLD  PRESSURE) 


A.  Reference  12 

Oi 

Average  Ni  (Drops/m3) 

AD  “  0 .5m* 

Range  of  Ni 

0.5 

18000 

1.0 

7500 

6900  -  2900 

1.5 

3200 

6900  -  1200 

2.0 

1400 

2300  -  650 

2.5 

600 

790  -  220 

3.0 

250 

460  -  ISO 

3.5 

100 

240  -  12 

lt.0 

50 

81  -  12 

4.5 

19 

5. 

a 

B.  Reference  13-  •  [Flour  Pellet  Sampling  Method  (33.1  tn/hr  fall).  Two  extreme  samples.] 


No.  Drops/m3 


Diameter  (mm) 


Sample  2 


Sample  A 


0.5  -  1.0 

5910 

4400 

1.0  -  1.4 

>279 

891 

1.4  -  2.0 

1247 

883 

2.0  -  2.4 

519 

407 

2.4  -  2.8 

kU 

419 

2.8  -  3-3 

V  7 

Vb 

3.3  -  4.0 

127 

139 

4.0  -  4.7 

44 

71 

4.7  ♦ 

18 

40 

N  -  9938. 


7566. 


TABLE  9.  SANDIA  FACILITY  DROP  DISTRIBUTION 
(TERMINAL  VELOCITY  NOT  ATTAINED)  FOR  RATE  ABOUT  5.5  IN/HR 
(OROPS  PHOTOGRAPHED  AT  OIL  INTERFACE) 


01  (m*) 

Ni  (Drops/ft3) 

Weiqht  Grams  (I0'3) 

0.5 

48 

3.23 

1.0 

90.7 

47.6 

1.5 

70.6 

125. 

2.0 

14.5 

60.8 

2.5 

0.7 

5-9 

30 

0.7 

10.2 

3-5 

0.7 

16.4 

N  -  226/ft3 

FORCE -PRESSURE -TIME  PROFILE  OF  IMPACTING  WATER  PROPS 

From  a  practical  aspect,  knowledge  of  the  dynamic  Interaction  of  a  water  drop  with  a 
finite  solid  surface  system  at  any  impact  velocity  appears  Incomplete,  both  analytically  and 
experimentally.  The  situation  appears  worse  when  the  "surface"  Is  a  system  which  has  signifi¬ 
cant  gross  motion  during  the  interaction.  Hydrodynamic  codes  are  being  employed  for  hypersonic 
collisions  of  solid  targets  with  solid  and  liquid  particle  project i les ^ ^  and  experiments  with 
piezoelectric  transducers  for  determining  force-time  from  supersonic  impacts  in  heavy  rainfall 
are  cont inui ng. ^ Several  complications  arise  In  the  general  question  of  the  specific 
force-pressure-time  profile  between  drop  and  system: 

a)  The  drop  is  not  always  spherical  in  a  rainfield. 

b)  The  bow  shock  preceding  a  supersonic  system  offers  some  degree  of  acceleration 
and  deformation  of  the  drop  before  system  collision. 

c)  Complicated  by  geometry  and  system  and  drop  compressibility,  a  pattern  of  pres¬ 
sure  and  tension  waves  traverse  the  drop  as  partially  illustrated  In  Figure  8. 


FIGURE  8.  QUALITATIVE  PROFILES  OF  INITIAL  WAVES  IN  LIQUID  DROP 

INELASTIC  COLLISION 


,,,,,, r-r,  yrrrrrr 

IVI  ]V|  I V, 


C=  COMPRESSION  WAVE  VELOCITY 

V|-  IMPACT  VELOCITY 

V=  PARTICLE  VELOCITY  BEHIND  SHOCK 
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At  first,  as  a  rigid  mass  at  velocity  Vj  strikes  a  stationary  drop,  a  shock  wave  originates  at 
the  interface,  moves  with  velocity  C  into  the  drop,  accelerating  the  water  traversed  to  some 
average  speed  V.  As  the  shock  propagates,  reflected  tension  release  waves  at  speeds  generally 
less  than  C  are  reflected  from  the  free  boundary  and  penetrate  the  shocked  fluid.  A  rarefac¬ 
tion  shock  cannot  exist  in  a  liquid.  A  pressure  gradient  Induces  radial  flow  while  water  above 
the  shock  front  remains  at  zero  velocity.  Some  spall inq  at  the  back  of  the  drop  might  be 
expected  as  the  reverse  tension  wave  is  initiated.  It  appears  that  the  drop  encounter  is  es¬ 
sentially  Inelastic  (or  plastic)  as  the  mass  of  the  drop  >s  splashed  in  a  direction  along  the 
surface  such  that  for  a  normal  encounter  as  illustrated,  the  momentum  imparted  to  the  system  « 
mVj.  Mo  experimental  references  cite  a  significant  rebound  of  the  drop  mass  on  a  rioid  sur¬ 
face.  Photographs  of  a  2  mm.  diameter  drop  on  a  hard  smooth  magnesium  alloy  at  1000  ft/sec 
revealed  negligible  spalling.  Correlation  of  energy  required  to  statically  deform  a 

honeycomb  crush  type  P0  fuze  with  raindrop  energy  available  from  high  velocity  sled  tests  at 
Holloman  Air  Force  Base  indicated  a  plastic  type  col  I  Is  ion. ^ 

d)  The  net  axial  force  at  any  time  is  the  Inteqral  of  pressure  over  the  effective 
contact  area;  but  pressure  appears  as  some  function  of  radial  distance  alonq  the  contact  area. 
Increasing  during  the  earlier  stages,  and  decreasing  durinq  the  later  stages  of  the  encounter. 

Some  liberty  will  be  used  In  interpreting  four  chosen  experimental-analytical 
models  to  arrive  at  the  magnitude  of  pressure  expected  on  PD  fuzes. 

Model  1.  Incompressible  Drop- I ncremental  Momentum  Transfer 
An  unrealistic  model  of  the  encounter  takes  no  consideration  of  shock  wave  propaga¬ 
tion  into  the  drop,  with  the  drop  impartlnq  forward  momentian  to  the  system  in  incremental  mass 
slices  dm  as  the  system  surface  with  velocity  V]  sweeps  over  the  drop  (Figure  9).  If  the  sys¬ 
tem  is  very  massive  and  rigid; 


FIGURE  9-  INCREMENTAL  MASS  MODEL 


dm  -irRoV,r(0-r)dt 

Force  F  ■  V|  dm/dt  «TR0V| 3t(D-Vj t) 
for  0  <  t  <  D/Vj 


(12) 

(13) 
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where  m  -  Drop  mass 


r  -  Distance  into  drop 

RQ  *  Density  of  uncompressed  water 

t  -  Time 

The  drag  coefficient  force  expression  with  the  dynamic  pressure  is: 

F  -  ed  *0  V,2  A/2  (H 

This  is  identical  to  expression  (13)  with  drag  coefficient  Cj  ■  2  for  plastic  encounter,  and 
the  contact  area  interface  A  »  »V|  t(D-Vjt).  Such  an  encounter  implies  a  constant  pressure 
R0 V | 2  throughout  the  collision  and  a  parabolic  profile  of  the  contact  force  (Figure  10). 
Representative  values  are  In  Table  10. 

FIGURE  10.  INCOMPRESSIBLE  DROP  FORCE-TIME  PROFILE 

F 


Peak  Force  «  *Ro(VjD)2 
V 


jV  dt  •  mV | 


-T  =  D/V,— *| 
m  =  DROP  MASS 


TABLE  10.  REPRESENTATIVE  VALUES  IN  INCOMPRESSIBLE  DROP  MODEL 


)  (rm) 

V,  (Ft/Sec) 

Pressure  (KSI) 

Peak  (Lb) 

T  (Microsec] 

2000 

54. 

262. 

3.28 

2 

3000 

121. 

590. 

2.18 

4000 

216. 

1050. 

1.64 

2000 

54. 

590. 

4.92 

3 

2000 

121  . 

1330. 

3.28 

4000 

216. 

2360. 

2.46 

2000 

54. 

2360. 

9.84 

6 

3000 

121. 

5310. 

6.56 

4000 

216. 

9450. 

4.92 

Mode)  2.  Col  I iding-Compres|ible  Cylinders,  Laterally  Constrained 
Though  neglecting  the  geometry  of  drop  and  system,  a  more  realistic  approach  for 
interpreting  initial  stages  of  contact  is  the  usual  column  of  height  H  of  compressible  water, 
constrained  laterally  to  avoid  pressure  release  effect,  struck  by  a  similar  infinitely  wide 
elastic  solid  surface  system  with  the  generation  of  planar  pressure  waves.  Fioure  II  shows  a 
cross  section  cut  out  of  an  infinite  width  of  the  process. 


FI  PURE  II.  INFINITE  PLANE  COMPRESSIBLE  WATER-SOLID  SYSTEM 

PRESSURE  BALANCE: 


WATER 

Poco 

C 


1. _ L 

/  /////  /  a / / 


(V^V) 


///////  /  / 


i 


INTERFACE 

/ 

7 


SYSTEM 
Pi  C! 


Vi 


P  =  PiCt  (Vi— VI  =  PqCV 


P  = 


PqCVi 

i+£o£' 
Pi  Cl 


PqCV! 


(WATERHAMMER) 


where 


(V, -V) 

zo 

zl 

Z 


-  Acoustic  velocity  In  water  (A800  ft/sec  at  20 ” C ) 

-  Acoustic  velocity  in  solid  system 

-  Shock  pressure  wave  velocity  in  water 

-  Height  of  water  column  initially  at  rest 

-  Density  of  uncompressed  water 

-  Density  of  solid  system 

-  Impact  velocity  of  solid  system 

-  Particle  velocity  in  water,  or  interface  velocity  durina 

-  Particle  velocity  in  system  surface 

-  Impedance  of  water  ■  R0CQ 

-  Impedance  of  solid  system  »  R|C| 


impact 
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For  pressure  balance  at  interface: 


P  »  R.C,(V,-V)  -  R  CV  -  z  vl 

'  '  r*  z/zT 


(15) 


in  Z,  the  wave  speed  cor  respond i nq  to  the  pressure  rise  created  by  the  wave  must  be  used  and 
not  the  acoustic  speed  of  the  undisturbed  liquid.  An  empirical  expression  has  been  reported 
for  the  shock  wave  velocity  C  at  which  the  pressure  wave  propaqates  upstream  throuqh  the 
f.uid:<l8> 


k  -  2  for  0  <  V/Cn  <  1.2 


k  =  I .33  for  I  .2  <  V/CQ  <  2 


(16) 


where  V  is  the  particle  velocity  in  water  and  is  approxi ma t e ly  expressed  as: 

:  V|B 


V  =  V 


2|/Z0  ♦  k (V) /C0)  1 

♦  (2,/Z0)  ♦  2k(V,/C0)  J 


on  the  condition  that: 

k(V,/CQ)  (I  +  ts(V,/C0)  ]  <<  U  ♦  (Z,/ZQ)  +  2k(V,/C0))2 
A  more  approximate  expression  for  V  is: 


V  «  V 


— ! - 1 3 

*  VhJ 


5  V,  B  ' 


(17) 


08) 


09) 


Thus,  the  pressure  becomes: 

P  =  R0(C0V  ♦  k V2)  (20) 

The  impedance  ratios  for  B'  for  aluminum,  steel-water  interfaces,  with  B  and  C  (ft/sec)  for 
various  impact  velocities  V|  are  listed  in  Table  II. 

TABLE  II.  IMPEDANCE  RATIOS,  PARTICLE  VELOCITIES,  SHOCK  VELOCITIES 
FOR  ALUMINUM,  STEEL-WATER  INTERFACES 


Aluminum-Water 

Steel- 

Water 

Zl'Zo 

9 

U 

26.7 

B' 

0 

•  9 

0.96 

V,  (Ft/Sec) 

B 

C 

B 

C 

100 

O 

k980 

.96 

5000 

1000 

.87* 

6500 

•  95 

6700 

2000 

.85 

8200 

•  91* 

8560 

J000 

.83 

9780 

■  93 

10400 

£888 

.81 
•  77 

1)280 

14000 

M 

12100 

15500 
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This  approach  indicates  shock  wave  velocity  much  higher  than  acoustic  velocity  in 


water  with  water  molecules  accelerated  close  to  the  initial  impinqement  velocity  V|.  The 
force  exerted  at  the  Interface  area  A  would  be  a  constant:  F  -  AP .  If  the  surface  system  is 
taken  to  be  incompressible,  Zj  qoes  to  infinity,  and  the  usual  waterhammer  expression  results: 

P  -  R0cv,  (21) 

as  the  pressure  existinq,  e.q.,  at  a  valve  when  suddenly  closed  to  water  movlnq  at  velocity  V| 
in  a  pipe.  If  the  water  is  of  heiqht  H,  this  pressure  exists  for  approximately  2H/C,  when  the 
reflected  pressure  release  wave  beqins  arrivino  at  the  valve,  resultinq  in  a  neqative  pressure 
contribution  there.  Dissipation  dampens  these  reflections. 

A  momentum  approach  results  in  the  same  waterhammer  yxrwr  .-n :  If  a  slab  of  liquid 
of  area  A,  thickness  x,  is  accelerated  to  a  velocity  V|  dur'^o  <*.<  me  7  it  takes  a  stress 
wave  at  speed  C  to  move  over  x,  then  the  force  on  the  slat.  «  ■-  r  -  R0ACVj.  The  pressure 
is  RQCV|.  Deviation  from  reality  of  the  compressible  cylinrh,  -  waachs 

1.  The  Surface  system  is  a  finite  solid  acte.  on  over  j  small  area,  with  corres¬ 
ponding  smaller  system  particle  velocity  (Uj  -  V)  than  is  'mi.iied. 

2.  The  pressure  qradient  between  the  inside  where  waterhammer  pressure  exists  and 
the  free  boundary  causes  radial  flow  at  velocity  Vr.  Initially,  reqardless  of  qeometry,  the 
radial  flow  stagnation  pressure  (R0Vr2/2)should  qrossly  equal  the  waterhammer  pressure 
(R0CV|)  or  Vr  » '^2CV|I.  Vr  becomes  several  times  the  impact  velocity  but  only  for  a  small 
portion  of  the  drop  mass  flow.  For  a  flat  ended,  free  boundary  liquid  cylinder  of  radius  R, 
the  impact  interface  area  over  which  the  waterhammer  pressure  acts  is  released  in  about  the 
time  R/Cc  after  impact.  The  drop  is  not  a  flat-ended  cylinder,  but  offers  chanqinq  contact 
area  as  collision  proqresses. 

Model  3-  Hydrodynamic  Incompressible  Drop  and  Rigid  Solid  Surface 

This  treatment ^ '  9)  concluded  that  maximum  pressure  occurs  In  an  outer  rlnq  of  the 
circle  of  contact  (Figure  12),  with  pressure  increasing  with  the  radial  distance  along  the 
plane  of  contact  up  to  this  rlnq.  When  the  ring  of  contact  is  about  O.tR,  this  rinq  of  maximum 
pressure  vanishes  and  maximum  pressure,  now  the  dynamic  pressure  R0Vj2/2,  exists  at  the  center 
with  pressure  decreasing  across  the  contact  area  from  the  center  to  the  periphery  of  the  circle 
of  contact  (Table  12). 
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FIGURE  12.  DROP-SURFACE  INTERFACE  TABLE  12. 


ity  of  a  spherical  drop  in  that  the  waves  are  not  all  started  simultaneously,  but  employs  some 
planar  wave  assumpt ions . (20,21 ,22)  |n  a  reoion  traversed  by  the  compression  wave,oiV|  is 
defined  as  the  average  axial  velocity  acquired  by  the  water  molecules,  of,  thus,  is  a  measure 
of  what  fraction  of  impact  velocity  V|  is  imparted  to  the  water  molecules  on  the  averaqe  from 
the  compression  wave  traveling  through  the  drop.  (of *  0  would  imply  no  compression  wave; 
of*  I,  that  the  wave  brings  water  molecules  up  to  the  striking  surface  velocity.)  oc  I s  said  to 
be  governed  mainly  by  divergence  of  the  wave  as  it  spreads  through  the  drop.  As  collision 
continues,  shock  waves  are  initiated  at  contact  points  of  the  sphere  and  the  X  plane  (upper) 
boundary  of  radial  flow  (see  Figure  13).  Maximum  pressure  may  be  expected  at  reqion  Q  where 
water  is  beinq  accelerated  both  by  the  shock  wave  arrival  and  by  the  effective  Ov  plane)  solid 
surface  front.  As  velocity  V)  increases,  wave  divergence  decreases  and  of  approaches  unity. 


FIGURE  13.  CONFIGURATION  (ENGEL)  LEADING  TO  TIME  TO  MAXIMUM  PRESSURE  AT  OUTER 


RING 


CT  (11) 
2R 

PRESSURE 
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With  0.9  taken  arbitrarily  for  o( ,  and  C  from  equation  lb  for  a  I  uni num-wa ter ,  values  for  this 


approach  are  listed  in  Table  13. 

TABLE  13.  TIME  T  TO  MAXIMUM  RING  PRESSURE  (tf  -  0.9,  ALUMINUM-WATER) 


D (mm)  1 

Vj (Ft/Sec! 

1  ')(0eq) 

T(Microsec) 

C  (Ft/Sec) 

2 

1000 

.9 

.03 

6900 

3000 

1 .8 

.09 

9780 

5 

1000 

.9 

.08 

6500 

3000 

1.8 

.10 

9780 

However,  it  has  been  reported  that  hiqh  speed  lateral  flow  bepins  when  0  -  17“ 
independent  of  impact  velocity;  and  that  maximum  impact  force  is  approximately  proportional  to 
impact  velocity,  implyina  constant  effective  Impact  area. 

It  is  stated  that  the  (averaqe)  maximum  pressure  created  durino  the  compression  wave 

period  is: 

P  -OfR0CV,/2  (23) 

This  expression  relates  to  the  time  aver aqed  impact  pressure  over  the  entire  collision  rather 
than  to  the  maximum  pressure  which  occurs  at  the  edaes  of  the  contact  area,  and  which  decays 
rapidly  because  of  radial  flow.  Values  for  o(  were  determined  from  pit-depth  equation  which 
involved  the  time  for  the  entire  process.  Some  inferences  about  were  derived  by  considerlnq 
that  particle  velocity  imparted  to  the  drop  shouH  be  some  fraction  e<  (to  maintain  contact  at 
the  interface)  of  the  velocity  chanqe  in  the  solid  system.  The  usual  constrained  cylinder 
intersects  the  periphery  of  the  drop  (fiqure  14), 


FIGURE  14.  CONFIGURATION  FOR  «  DERIVATION 

V  <  (Vj-W) 

V  °f(Vj  -W)  (24) 

•  >  _  ,  (I-*)  V|  i  Vj-V-** V  ar  V | - V  (for  riqld  surface) 

I  I  w  ‘ 

*  v  <y«v, 

where  V  -  Particle  velocity  In  compressed  water 

Vj  -  Initial  Impact  velocity  of  solid  on  stationary  drop 
W  -  Particle  velocity  In  compressed  solid 
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V*»v- 


I 


{ I -o*) V |  becomes  the  velocity  at  which  the  surface  moves  through  the  drop-  This  appears  to 
effectively  lengthen  the  collision  duration  by  a  factor  l/(l-*t)  lonqer  than  D/V j  .  In  an 

experiment  with  a  5.7  mm  diameter  drop  at  Its  V|  -  26.9  ft/sec  terminal  velocity,  the  value  of 
was  determined  by  observing  the  radial  flow  velocity  and  using  the  Bernoulli  relation 
Vr  »  >/«<  CVj  to  determine  «f  as  0.4.  A  barium  titanate  disk  coated  with  silver,  lacquered  and 
cemented  with  polystrene  adhesive  to  a  metal  base  showed  the  force  for  a  40  foot  fall  of  this 
drop  rose  very  rapidly  and  underwent  a  rapid  decay  to  zero  in  about  one  millisecond  (figure 
15).  The  interaction  time  appeared  extended  from  0.7  ms  ( D/V | )  to  about  1.2  ms. 

FIGURE  15.  fORCE  PROFILE  FROM  LOW  VELOCITY  DROP 


It  appears  difficult  to  set  up  force  and  pressure  time  curves  for  any  of  the  last 
three  models.  Calculated  comparison  pressure  magnitudes  for  the  four  models  are  listed  In 
Table  14. 


TABLE  14.  COMPARISON  OF  ESTIMATED  INTERFACE  PRESSURE  (KSl) 
OF  VARIOUS  MODELS  FOR  5  MM  DIAMETER  WATER  DROP 


Aluminum-Water 

Steel 

-Water 

Impact  Velocity  V| (Ft/Sec) 

1500 

3000 

1500 

3000 

1. 

R0VjJ  (Momentum  Transfer) 

30 

121 

30 

121 

2. a. 

RcCV|  (Modified  Waterhammer) 

1  ♦  «o5_ 

R 1 C  i 

128 

324 

145 

386 

b. 

RqCV|  (Waterhammer) 

149 

394 

154 

418 

3. 

3R(jV)2  (Hydrodynamic) 

91 

364 

91 

364 

4. 

RqCV,/2  (Engel) 

75 

197 

77 

209 

Water  shock  velocity  C  calculated  from  (15) .  of(Engel)  taken  as  I.  Momentum  trans¬ 
fer  model  pressure  constant  for  entire  collision.  Other  pressures  are  pertinent  for  "Initial" 
stages  of  collision. 
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SIMULATION  OF  SPRING-MASS  PO  AT  HIGH  VELOCITY  IN  RAIN 


For  an  analytical  description,  the  system  impacted  appears  as  an  involved  mechanical 

system.  Three  common  rain  PO  desensitizers  (Fiqure  16)  are  Illustrated:  (I)  Crusher  wl th  a 

honeycomb  fixture  supportinq  the  PO  for  absorblnq  the  kinetic  enerqy  of  drops:  (2)  Buffer  wi th 

some  effective  sprinq  mass  system  which  buffers  momentum  onto  the  shell,  continuously  restor- 

(2l|) 

inq  Itself  to  Initial  position;  (3)  S I  I cer  employinq  some  recessed  cavity  with,  e.q., 
crossbars  of  sufficient  strenqth  to  shatter  the  drops  Into  smaller  size  and  simultaneously 
buffer  some  momentum  onto  the  shell.  Soft  tarqet  sensitivity  and  erosion  dictate  limits  on 
Such  des i qns  . 


FIGURE  16.  PO  RAIN  DESENSITIZING  DESIGNS 


2.  BUFFER  (ELASTIC)  3.  SLICER  (CROSSBARS) 


Simulation  Procedure  for  Buffer  System 


In  an  elementary  approach,  the  conservation  of  momentum  mv|  (m  -  drop  mass,  vj  • 
impingement  velocity)  can  be  used  with  certainty,  and  the  initial  enerqy  with  reservation 
because  of  Its  uncertain  distribution  between  consequent  gross  system  motion  and  heat  dissipa¬ 
tion;  but  neither  is  of  immediate  assistance  in  the  time  domain.  The  response  of  a  sprlnq- 
mass  PD  (Fiqure  17)  in  a  proposed  artillery  fuze  to  hiqh  velocity  fliqht  In  heavy  rainfall 
was  simulated  by  a  straiqht-forward  approach  employinq  conservation  of  momentum  and  ordinary 
harmonic  motion  for  the  time  between  raindrop  impacts.  A  sprinq,  loosely  positioned  in  a 
confining  channel  In  the  steel  PD  nose,  places  the  system  under  a  40  pound  nominal  compres- 
slonal  load,  and  a  further  compression  of  0.26  Inches  is  needed  to  bring  the  firing  pins  In  a 
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position  to  begin  Initiating  sensitive  M55  stab  detonators.  If  the  mass  and  spring  were 
attached  and  free,  the  natural  period  +  would  be  14.5  ms.  For  velocities  of 

interest  in  the  range  1000  to  3000  ft/sec,  the  expected  force  contact  with  any  drop  Is  In  the 
order  of  several  microseconds;  and,  therefore,  the  force-time  curve  (unknown)  was  avoided  and 
the  momentum  I  transferred  or  area  under  the  force-time  curve  was  employed. 

I  «  Cf  dt  “  mv  ■  HV  (25) 


FIGURE  17.  MASS-SPRING  P0  DESIGN 
SPUR 


System  Parameters; 

A  -  Area  of  flat  front  end  -  .001364  ft2  (1/2"  diameter) 

D  -  Drop  diameter 

K  -  Spring  constant  (40  Ib/in  nominal) 
m  -  Mass  of  raindrop 

m,  -  Hass  of  spring  »  .00054  slug  (0.28  oz) 

M  -  Mass  of  steel  PD  nose  ■  .00235  slug  (1.21  oz) 

S  -  Compress  Iona)  static  load  (40  lb  nominal) 
v  -  Velocity  of  shell  relative  to  ground  (ft/sec) 

V  -  Velocity  of  PD  nose  relative  to  shell  body  (ft/sec)  ■  x 
w  -  Natural  freguency  of  free  system  •  435  rad/sec 
x  -  Position  of  PD  relative  to  shell  body 
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it 

With  harmonic  motion  between  impacts  and  neqlect  of  all  friction,  especially  In  the 


0  ring  area,  the  equation  of  motion  and  conditions  are  simply: 

x  (0)  -  Xj  , 

(M  +  m^/3)  x  -  -Kx  x  (0)  »  V  (0)  •  Vj|  Constraint:  I "  <  x  1.26"  (26) 


with  position  x  and  velocity 

V  of  the  PD  nose : 

x  (t)  *  V; |  sin  wt 

♦  Xj 1  cos  wt 

(27) 

w 

V  (t)  «  Vj |  cos  wt 

-  w  x  j  j  sin  wt 

(28) 

xmax  *  xi 

♦  (Vn/w  Xj,)2' 

^  Valid  for  +  V . | 

(29) 

TVELO  ■=  l  tan 

(Vj,/w  Xj | )  J 

1  (Motion  toward  detonators) 

(30) 

w 


Xj  -  PD  position  just  BEFORE  encounter  with  ith  drop. 

Xj .  -  PD  position  just  AFTER  encounter  with  ith  drop  =  «j  ♦  VjDj/v. 

Vj  -  PD  velocity  just  BEFORE  encounter  with  ith  drop. 

Vjj  -  PD  velocity  just  AFTER  encounter  with  ith  drop  *  Vj  +  mjV/M. 

xmax  '  MaKinu,r  possible  excursion  after  any  impact. 

TVELO  -  Time  for  xmax  to  occur. 

Calculated  values  for  the  PD  design  for  sirioie  drop  encounters  are  listed  in  Table  15. 


TABLE  15.  PARAMETERS  FOR  SI  HOLE  DROP  ENCOUNTER  FOR  SYSTEM  (S«A0  lb;  K-AO  Ib/in) 


Drop  Diameter 
D(mm) 

Shell  Velocity 
v(ft/sec) 

1 mpu 1 se 

1 (o? -sec) 

PD  Ve loc i ty 
V(ft/sec) 

xmax^nr’) 

TVELO (ms) 

2 

2000 

0.0071 

0.19 

•v  0 

_ 

3000 

0.0122 

0.33 

n.  0 

- 

3 

2000 

0.027 

0.73 

0.005 

0.05 

3000 

0.0AA 

1  .18 

0.013 

0.07 

5 

2000 

0.138 

3.68 

0.13 

0.23 

3000 

0.212 

563 

0.30 

0.35 

6 

2000 

0.237 

6.31 

0.38 

0.39 

_  . 

3000 

0.368 

9.78 

0.91 

0.60 

7 

2000 

0.386 

10.3 

'  1.0 

0.63 

3000 

0.586 

15.6 

2.25 

0.93 

8 

A  000 

1.16 

31 .0 

8.0 

1.62 
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It  was  assumed  drops  were  uniformly  distributed  in  space,  and  the  time  between  col¬ 
lisions  (TBCOl)  was  determined  by  the  volume  swept  out  by  the  PD  area  A,  with  a  shell  velocity 
v,  in  rain  density  N.  TBCOL  *  l/(NAv).  Two  distributions  were  used:  I)  the  cited  Canal  Zone 


114  mm/hr  (4.5  In/hr)  equivalent  to  a  cloudburst  (liquid  content  4.5  g/m3)  ;  and  2)  Holloman 
"deluge"  ralnfield  of  about  28  in/hr  accumulation  rate  (liquid  content  25  q/m3)  from  a  flour 
pellet  measurement  with  the  number  of  drops  at  the  tarqe  diameter  end  arbitrarily  Increased 
to  exagqerate  (liquid  content  33  g/m3)  this  deluge  rainfield  for  a  more  pessimistic  response 
Each  successive  drop  was  chosen  randomly  from  the  drop  size  distributions  (Table  16). 

TABLE  16.  RAIN  FIELDS  USED  IN  SIMULATION 


Canal 

Zone  (Table  4;  11 k  mm/hr) 

Exaggerated  Holloman 

Del uqe 

Nj/m3 

l  Total  Drops 

D(mm) 

Nj/m3 

$  Total  Drops 

D(nin) 

141 

9-7 

•  7 

4400 

58 

•  73 

373 

25-8 

1 

891 

12 

1  .26 

281 

19-3 

1.3 

883 

12 

1.77 

204 

14.0 

1.6 

407 

5-1 

2.31 

147 

10.1 

1.9 

419 

5-5 

2.7 

128 

8.8 

2.2 

316 

4.2 

3.42 

(3.1) 

80 

5-5 

2.5 

139 

1.8 

4.32 

(3-7) 

46 

3.2 

2.8 

71 

.9 

5.04 

(4.4) 

25 

1.7 

3.1 

40 

.5 

6.2 

(4.7+) 

14 

■  96 

3-4 

N  -  7566 

100 

7-5 

.52 

3-7 

3-5 

.24 

4 

(Values  in 

parenthesis  for  Holloman  were 

1.4 

.1 

4.3 

' averaqe  ' 

diameter  values  for 

actual  4th 

0.6 

.0' 

4.6 

sample  In 

Table  8.) 

0.3 

.02 

4.9 

0.3 

.02 

6 

N  -  1452.6 

100 

In  the  simulation  procedure,  final  conditions  of  position  and  velocity  of  the  PD  just 
before  another  drop  impacted,  along  with  an  impulsive  velocity  from  the  next  drop,  became  the 
initial  conditions  for  harmonic  motion  in  the  next  TBCOL.  Whenever  Xj  £  0.2  mu,  printout  occur¬ 
red  of  Xj  and  Vj  at  each  T.  (Tj  *  TBC0L/5) .  All  drop  sizes  were  printed.  Ordinary  shell 
velocity  trajectories  were  employed,  and  also  velocities  at  4000  ft/sec.  For  the  realistic 
case  of  3000  ft/sec,  simulation  indicated  neqliqible  PD  motion  for  the  Canal  Zone  distribution, 
and  excursions  up  to  about  one  mm  for  the  exaggerated  Holloman  deluge  distribution.  Figure  18a 


is  a  typ  cal  computer  plot  of  PD  displacement  toward  the  detonators  (about  7  mm  needed)  versus 
distance  travel  (2000  ft)  In  the  ralnfleld  and  shows  the  recovery  feature  of  the  nose.  Subse~ 
quent  to  the  simulation,  prototype  fuzes  were  built  and  tested  (December  1969)  at  Holloman  AFB 
in  the  deluge  rain  (=28  in/hr)  with  an  entrance  rocket  sled  velocity  of  about  2750  ft/sec  Into 
the  2000  ft  of  rain.  Inspection  of  fuzes  revealed  negligible  maximum  motion  (small  fraction  of 
millimeter)  of  the  PD  nose. 

m  FIGURE  18.  SIMULATED  P0  MOTION  (PD  DISPLACEMENT  MM  VERSUS  SHELL  TRAVEL  FT) 
f  a.  Exaggerated  Holloman  Deluge  (liquid  content  33  g/m3) 


7I  DET  LINE -_7_rri 


v  «  3000  ft/sec 
TBC0L  -  1.14  ms 
DROPS  IMPACTING  -  580 


b.  Truncated,  Resonant  Holloman  Deluge  (liquid  content  •  2870  g/m3) 


DET  LINE”  7  mm 


HOLLOMAN  (TRUNCATED)  —  Drops  3.4  to  6.2  n»n 
v  »  2000  ft/sec 
TBC0L  -  0.2  ms 

DROPS  IMPACTING  -  500  ^ 

c.  Truncated,  Resonant  Holloman  Deluge  (liguid  content  •  1910  g/m3) 


v  •  3000  ft/sec 
TBC0L  -  0.2  ms 
DROPS  IMPACTING  -  400 


- — -  240  feet 

Period  -  9.1  ms 
(Natural  Period  «  14.6  ms) 

"Reflection  of  neqative  peaks  gives  semblance 
of  2nd  harmonic,  heavily  damped." 


Disregarding  the  surge  wave  problem,  the  above  approach  was  employed  for  a  "trun¬ 


cated"  Holloman  field  by  disregarding  small  size  drops,  then  increasing  this  distribution  for 
a  "resonant"  condition  of  drops  Impacting  the  PD  every  0.2  ms.  This  effectively  increased  the 
liquid  content  of  the  Holloman  deluge  by  about  100.  At  this  intensity,  the  total  water  content 
is  still  only  a  fraction  of  one  percent  of  the  total  volume.  Table  17  lists  these  values, 
figures  18b  and  c  are  typical  computer  outputs  for  these  bizarre  type  rainfields. 


TABLE  17.  ARBITRARY  LARGE  DROP.  HIGH  DENSITY  RAINFIELDS 


Truncated  Holloman 

"Resonant"  Holloman  (T8C0L  ■  0.2  mil llsec) 

D(rnti) 

Nj/m3 

N,/m3(v-2000) 

N|/m3(v"3000  ft/sec) 

3 . 

193 

28,230 

18,000 

k. 32 

139 

20,330 

13,550 

5-04 

71 

10,360 

6.900 

6.2 

(o 

5.830 

3.880 

N 

443 

64,750 

42,330 

Liquid  Content  (q/m3) 

197 

2870. 

1910. 

Criticisms  of  this  Hethod 

I.  The  assumptions  of  simple  harmonic  motion  between  impacts,  and  of  conservation 


of  enerqy  In  the  PD  mass  from  the  plastic  drop  encounters  until  the  PD  mass  strikes  Its  mech¬ 
anical  stop  at  x  •  one  inch,  at  which  all  acquired  momentum  and  enerqy  is  taken  as  transferred 
to  the  shell  itself,  take  no  account  of  wave  propagation  effects  alonq  the  spring.  There  is 
room  for  a  variety  of  alternate  approaches.  The  approach  of  stress  wave  propagation  in  the 
metal  nose,  the  boundary  conditions  on  the  spring  ends  and  nose,  with  momentum  buffering  onto 
the  shell  system,  the  standard  coupled  differential  equations  for  the  spring  with  associated 
masses  and  constants,  implies  Its  own  inherent  assumptions  and  offers  an  Interesting  approach 
of  much  greater  complexity  for  computer  solution.  For  a  spring  constrained  at  both  ends,  surge 
waves ^5)  occur  If  force  loadinq  occurs  at  frequencies  fn  (Table  18)  for  this  design. 

fp  «  nd  \J  flg  n  “  1,2,3  ...  (31) 

2*fDzN  V  2w 

d  -  Diameter  wire  m  .065"  N  -  Number  active  colls  *  20 

D  -  Diameter  of  coll  «  .285"  w  -  Weight  density  ■  0.286  Ib/ln3 


TABLE  18.  SURGE  WAVE  FREOUFNCIIS 


n 

fn(HZ> 

Period  (ms) 

1 

160 

6.3 

2 

320 

3.1 

3 

480 

2.1 

L 

( 26) 

A  reference'  cites  the  time  for  a  compress ionat  wave  to  travel  down  and  back  alonq 
a  sprinq  compressed  or  extended  to  a  lenoth  L  as  2  m^/K  '*  2.1  ms.  The  propagation  velocity 
for  the  design  is  L  \J~ K/ms  -  150  ft/sec  for  the  L  ■  1.88  inch  40  ib  load  compression  lenoth. 

2.  The  value  of  the  aerodynamic  force  from  air  flow  was  not  known.  This  force  makes 
the  PD  more  rain  sensitive,  subtracts  from  the  40  Ib  static  load,  and  shifts  the  sprinq  force- 
compression  curve  to  a  lower  parallel  line  requirinq  direct  use  of  equations  of  notion  (Flqure 
19).  (The  dynamic  air  pressure  term  RQv2/2  at  3000  ft/sec  for  this  nose  is  about  13  Ib.) 


FIGURE  19.  FORCES  OH  PD  SYSTEM 


T  *- 


jnrrdinnr(wl  . — *  5  ♦  * *■ 


T  -  Aerodynamic  Force 


(h  ♦  ms/3)  x  -  T  -  (S  ♦  Kx) 
x  (0)  -  0  0  <  x  <  0.26" 
x  (0)  »  Vj  | 


(32) 


Sensitivity  Curve  for  the  PD  Fuze 

The  enerqy  required  to  compress  the  spring  another  0.26  inch  (1.08  ft-lb)  [sprinq 
enters  sliqhtly  non-linear  region),  and  to  initiate  the  M55  detonators  in  this  design 
(  ’  .05  ft-lb),  requires  an  impulse  of  I. 16  oz-sec,  or  impulsive  PD  velocity  of  3'  ft/sec.  To 
obtain  an  indicative  sensitivity  curve  (Fiqure  20),  based  on  rectanqular  pulses,  a  plot  of  con¬ 
stant  force  level  Fq  that  must  be  maintained  for  time  Tq  is  plotted  versus  impulse  I Q (”FoT0) 


delivered.  With  the  equation  of  motion  (M  ♦  ms/3)  x  *  FQ  -  (40  +  Kx) , 
tlon  are  satisfied  when  at  time  TDET,  the  solution  attained  the  values 
sec.  Aerodynamic  loading  will  shift  the  asymptotes. 


conditions  for  detona- 
x  *  6.6  inn  and  x  •  6  ft/ 
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FIGURE  20.  CONSTANT  FORCE  —  PO  SPRING  HASS  SENSITIVITY  CURVE 


EROSION  PHENOMENA 

In  the  1000  to  3000  ft/sec  velocity  ranqe  erosion,  rather  than  cratering,  Is  the 
dominant  form  of  damage  from  raindrops  ^?)  and  all  structural  materials  show  erosion  damage 
when  exposed  to  nominal  rain  for  sufficiently  lonq  time  at  velocities  In  excess  of  Mach  0.8. 

In  1961,  an  F106A  airplane  flown  throuqh  thunderstorms  from  15,000  to  40,000  ft  from  Mach  0.84 
to  I.63  sustained  severe  erosion  from  water  drops  and/or  tiny  Ice  crystals  at  rivet  heads  and 
leading  edge  of  winqs,  cockpit  frame  and  plastic  antenna  covering.  The  principal  causes  cited 
for  erosion  damage  are:  (I)  impact  force;  (2)  radial  flow;  and  (3)  cavitation  (Fiqure  21). 
However,  there  are  almost  as  many  mechanisms  of  erosion  failure  as  there  are  classes  of 
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materials.  Thin  resilient  coatings  such  as  neoprenes  end  polyurethanes  will  transmit  shear 
stress  to  the  substrate  causinq  failure  In  the  adhesive  bond.  Additional  impacts  cause  the 
coatinq  to  stretch,  deform  or  burst.  Plastic  type  materials  and  soft  metals  flow  plastically 
resulting  In  cratering  and  pitting.  Plastic  laminated  materials  fail  from  erosion  of  upper 
layers  of  the  fabric. 

FIGURE  21.  CITED  CAUSES  OF  SURFACE  EROSION 

1.  IMPACT  FORCE  2.  RADIAL  FLOW  3  CAVITATION 


impact  Force 

During  the  early  staqes  of  collision  of  a  drop  with  a  planar  solid  surface,  maximum 
pressure  appears  independent  of  drop  diameter  and  is  asserted  to  exist  in  a  rina  around  the 

central  point  of  collision.  At  supersonic  velocities,  local  pressure  of  about  0.1  to  I  million 

(29) 

psi  is  exerted  on  the  target.  Part  of  the  Impact  enerqy  from  liquid  Impact  radiates  into 
the  solid  by:  (I)  lonqltudinal  compression  wave;  and  (2)  transverse  or  shear  wave  (at  about 
half  the  lonqltudinal  velocity)  with  particle  motion  perpendicular  to  propagation  direction. 

If  the  waves  reflect  on  a  free  surface  of  the  solid,  each  incident  wave  generally  produces  two 
reflected  waves,  with  a  complicated  stress  pattern  resultinq.  Large  tensile  stresses  are  pro¬ 
duced  when  Intersecting  waves  are  In  opposition,  and  lead  to  fractures  especially  in  brittle 
materials,  rather  than  In  ductile  metals.  The  spalling  feature  occurs  when  an  initial  compres¬ 
sion  wave  is  being  reflected  In  the  opposite  surface  as  a  tension  wave.  Fracture  occurs  a 


little  way  in  from  the  rear  surface  where  the  algebraic  sum  of  reflected  tensile  stress  and 
compressive  stress  associated  with  the  tail  end  of  the  incident  wave  first  reaches  the  solid 


tensile  strength.  The  volume  of  material  between  fracture  and  rear  surface  still  has  high  par¬ 
ticle  velocity  in  the  outward  direction  and  may  fly  off  the  back  as  a  spall.  The  fracture  sur¬ 
face  created  acts  as  a  new  reflecting  surface  for  the  remainder  of  the  compression  wave,  and 
may  lead  to  secondary  fractures  in  toward  the  center  of  the  plate,  fractures  near  the  Impact 
surface  are  caused  by  tensile  stresses  accompanying  the  release  of  the  load.  For  angles  other 
than  normal  Incidence,  the  energy  of  the  incident  compression  wave  reflected  In  the  rear  sur¬ 
face  is  shared  between  a  wave  of  tension  and  a  shear  wave. 

Radial  Flow 

A  solid  sphere  inflicts  damage  by  exerting  localized  pressure;  whereas  a  liguld  drop 
can  cause  damage  by  its  radial  flow,  exerting  an  erosive  shear  stress  on  the  solid  surface 
because  of  the  liguid  viscosity,  and  also  at  the  base  of  any  protrusion  on  the  surface  with  a 
bending  moment  about  the  protrusion.  As  the  impacting  drop  spreads  out,  the  water  radial 
velocity  may  be  two  to  three  times  impact  velocity.  For  a  glass  plate  surface,  flow  velocity 
approaching  ten  times  the  impingement  velocity  for  short  times  after  Impact  has  been  found. 

Table  19  lists  some  results  of  single  Impact  studies  with  8  million  frames/sec  camera 
rate  of  the  radial  velocity  of  2  mm  water  drops  at  987  ft/sec  onto  Perspex. 

TABLE  19.  2  nin  WATER  DROPS  ON  PERSPEX  (PLASTIC) 


Impact  Velocity  (ft/sec) 

Time  from  Impact  (Microsec) 

Radial  Velocity  (ft/sec) 

0.8  -  2.9 

3056 

987 

2.9  -  3-9 

2650 

3.9  -  5.6 

2610 

5.6  -  12.9 

1820 

No  spalling  effect  was  observed  at  the  back  of  the  water  drop.  However,  at  low  velo¬ 
cities  and  at  1000  ft/sec,  there  is  evidence  at  the  last  stages  of  collision  when  the  drop  has 
been  flattened  Into  a  radially  flowing  disk,  that  a  spray  of  much  smaller  droplets  at  the 
periphery  appears  to  rebound  from  the  striking  surface. 


123 


1  • 


Hard  spheres  do  not  satisfactorl ly  simulate  water  impact.  Polytetraf luoroeth 'lene 


(PTFE)  and  modeling  clay  (plasticine)  most  closely  resemble  water  flow  characteristics  but  have 
scourinq  action.  Table  20  shows  some  results'^  '  with  solid  spheres  which  showed  some  perma¬ 
nent  set,  except  for  sapphire  which  had  complete  recovery. 

TABLE  20.  SOLID  SPHERE  COLLISION  ON  PERSPEX 


Sphere  Material 

Contact  Duration 

D(mm) 

(ft/sec) 

(ml crosec) 

(ft/sec) 

Po 1 yethy 1 ene 

1018 

52. 

115 

A 

Ny  Ion 

957 

6.7 

2k0 

k 

Cellulose  Acetate 

985 

5. 

250 

L 

Sapph i re 

1075 

8.5 

373 

2 

Cavi tat  ion 


Cavitation  is  the  formation  of  bubbles  in  a  liquid,  and  occurs  when  pressure  in  a 
liquid,  or  in  a  small  volume  in  a  liquid,  drops  below  the  liquid  vapor  pressure  at  the  tempera¬ 
ture  in  question.  When  pressure  on  the  liquid  is  raised,  or  when  bubbles  move  out  of  a  local 
low  to  a  hiqh  pressure  reqion,  the  bubble-cavities  collapse.  Collapse  of  cavities  produces 
damage  called  cavitation  erosion.  Ourinq  the  duration  of  each  impact,  some  bubbles  will  be 
collapsing  on  or  near  the  impact  surface.  The  bubble  collapse  lasts  only  a  few  microseconds 
and  the  collapse  speeds  are  estimated  to  be  qreater  than  impact  velocity.  Shocks  from  implo¬ 
sion  offer  a  small  contribution  to  damaqe ,  while  the  main  damaqe  is  caused  by  jet  impact  near 
the  end  of  collapse  of  the  spherical  bubbles  in  an  asymmetric  mode  such  that  high  velocity 
microjets  impinqe  on  the  surface.  It  appears  cavitation  microjet  diameter  Is  well  less  than 
one  mil  and  with  velocity  up  to  several  thousand  ft/sec,  and  probably  with  large  length  to 
diameter  ratio.  ^8,JI  ,32  •  33 ) 

There  are  theorles^^  for  the  origin  of  this  cavitation:  (I)  When  the  head  of  the 
drop  has  just  disappeared  into  radial  flow,  the  continuous  outward  flow  of  liquid  under  its  own 
momentum  produces  a  pressure  drop  at  the  center  of  the  spreadinq  liquid  disk.  If  pressure  at 
the  center  of  the  spreadinq  liquid  d i • k  falls  below  the  vapor  pressure,  cavitation  results. 

(2)  Alternating  waves  of  compression  and  tension  exist  In  the  head  of  the  liquid.  The  initial 
compression  wave  Is  reflected  from  the  free  I iqul d- to-ai r  surface  as  tension  or  negative 


pressure  wave  and  moves  down  to  the  surface.  Returninq  tension  wave  adds  alqebraical ly  to  the 
compression  wave  still  beinq  initiated  at  the  impact  surface,  and  as  the  impact  pressure  b*s 
been  decreasinq  steadily,  the  net  pressure  may  be  neqatlve.  Repeated  reflection  of  the  wave 
may  occur  before  the  bubble  qrows  to  sufficient  size. 


Erosion  Factors 

Generally,  erosion  varies  as  the  cosine  of  the  ancle  of  incidence.  Sliqhtly  higher 

(18) 

erosion  rate  was  found  at  ten  deorees  from  the  normal'  '  for  aluminum  at  Mach  1.2:  and  for 

neoprene  coatinqs,  most  erosion  occurred  at  impact  anoles  between  25  and  30  deqrees  from  the 
(31*) 

normal.  With  drops  of  similar  size,  flattened  droplets  caused  more  severe  erosion  than 

oval  shaped  ones.  If  erosion  E  is  defined: 


E  „  Volume  of  material  lost  per  unit  area  per  unit  time 
Volume  of  liquid  impinged  per  unit  area  per  unit  time 

an  empirical  formulation  for  E  (Expression  33)  shows  the  stronq  dependence  on  velocity.  It  is 

also  reported  that  at  velocities  up  to  about  1300  ft/s ec,  weiqht  loss  from  samples  In  rain  Is 

proportional  to  the  fifth  power  of  velocity. 

E  5=  V3 (V  -  Vc)/Vc‘*  for  I  <  V/Vc  <  3  (33) 

V  -  Velocity  of  impinqement 

V  -  Cut  off  velocity  for  erosion  (e.g.  *  390  ft/sec  for  0.66  mm  drops) 


Maximum  impact  stress  generally  is  a  function  of  material  properties  and  impact  velo¬ 
city.  Vet,  erosion  also  depends  on  drop  size  and  shape.  Though  contact  area  varies  as  D2  for 
a  spnere  and  as  D  for  a  cylinder  jet,  the  force  from  cylindrical  jets  of  water  on  barium  ti- 
tanate  transducer  showed  loads  rather  linear  with  velocity,  without  a  clear  dependence  on  jet 
diameter. 

Impact  stress  alone  is  not  a  sole  criterion  for  erosion.  Larger  drops  increase  the 
time  of  force,  and  materials  with  definite  yield  point  in  the  static  stress-strain  curve  have 
measurable  time  delay  associated  with  the  Initiation  of  plastic  deformation.  Strain  rate  for 
the  usual  longitudinal  impact  test  is  particle  velocity  divided  by  specimen  lenqth.  For  impact 
on  a  semi  -  inf  I ni te  body,  lenqth  could  be  related  to  the  depth  of  quasi-static  stress  fields. 
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smalt  for  smaller  drops,  resultinq  In  high  effective  strain  rate.  If  there  Is  a  strain  rate 
effect,  such  that  for  a  given  strain  the  stress  increases  with  strain  rate,  then  maintaining 
the  same  total  volume  of  Impinqing  water,  the  total  strain  and  relative  energy  transfer  for  a 
given  impact  velocity  would  be  less  with  smaller  drops. 

Aerodynamic  heatinq  can  reduce  the  strain-stress  curve,  resultinq  in  damage. 
Polyurethane  was  found  to  be  auout  one-sixth  as  stronq  at  250“F  than  at  Jk*F  If  the  im¬ 

pact  stress  does  not  exceed  the  coating  tear  strength,  these  temperature  induced  reductions  of 
impedance  may  be  beneficial. 


Erosion  Abatement 

Various  coatings  and  coverings  are  employed  for  rain  protection.  Nickel  plat- 
ing^®’^  has  extended  life  of  the  leading  edge  materials  by  about  40  times  over  neoprene.  A 
thickness  of  at  least  12  mil  nickel  on  laminates  extended  life  of  leadinq  edges  on  aircraft, 
helicopter  rotor  blades  and  jet  enqine  compressor  blades,  and  a  16  mil  thickness  has  been 
employed  on  high  strenqth-temperature  plastic  (polybenzimidazole  PBI).  An  explanation  for 
improved  performance  from  a  rough  nickel  surface  asserts  that  rough  (sandblasted)  surface 
breaks  the  drop  into  minute  particles  and  these  minute  particles  produce  less  radial  flow  and 
'ower  shear  stresses. 

Limited  data  is  available  on  rain  erosion  at  velocities  from  1000  to  3000  ft/sec. 
Cushioning  of  the  impact  stresses  only  occurs  once  a  critical  thickness  of  coating  is  reached 
of  several  millimeters.  As  the  impact  pressure  reduces,  the  velocity  of  radial  flow  is  also 
reduced  with  accompanying  shear  surface  stress.  No  known  elastomeric  coating  or  glass  rein¬ 
forced  plastic  can  withstand  erosive  action  of  rain  at  velocities  greater  than  3000  ft/sec  for 
more  than  a  few  seconds.  Ceramic  materials  resist  damage  to  <t000  ft/sec.  Epoxy-glass  nose- 
cones  with  radius  of  about  1.5  Inches  eroded  after  85  seconds  exposure  at  1150  ft/sec;  8  sec¬ 
onds  at  2200  ft/sec;  and  3  seconds  at  2600  ft/sec  (Holloman  AFB  6000  ft  field).  Unprotected 
ceramic  nose  cones  shattered  in  the  sled  tests.  Composite  materials  do  not  withstand  liquid 
impact  well  because  of  the  sharp  stress  peak  caused  by  the  initial  compressible  behavior  of  the 
material,  and  stress  waves  produced  lead  to  Interlaminar  failure.  A  possible  composite  system 
advantage  is  the  attenuation  or  dispersion  of  the  stress  waves  by  multiple  reflection  processes. 
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If  the  dynamic  stress-strain  curve  is  plotted  for  some  material  from  data,  the  slope 
yields  the  elastic  modulus  E,  and  the  elastic  wave  speed  C  is  then  (E/R)^'^  where  R  Is  the 
density.  Dynamic  impedance  Z  -  RC .  Elastomer  (several  polyurethanes,  neoprene)  coatinqs  have 
small  impedances,  and  in  effect  qlve  easily  under  raindrop  Impact.  These  materials  can  endure 
large,  reversible  amount  of  deformation  without  beinq  damaged .  Elastomeric  materials  resist 
erosion  (at  prescribed  temperature,  and  except  after  extreme  prolonged  exposure)  by  minimizing 
stresses  at  impact  area,  provided  the  velocities  are  sufficiently  below  the  level  needed  to 
fracture  the  elastomer.  The  material  must  recover  fast  enough,  retainlno  its  elastic  proper¬ 
ties.  A  random,  hiqh  strain  rate  loadinq  exists  from  rain  impact.  The  dynamic  response  of  a 
polyurethane  at  a  strain  rate  of  about  1000  sec"'  showed  stresses  three  times  higher  (at  com¬ 
parable  strains)  than  the  stresses  obtained  statical  ly. 

Experimental  studies  of  polyurethanes  have  indicated  that  low  modulus  with  qood  ten- 

11,  | ) 

sile  strenqth  and  elonqation  qlves  the  most  erosion  resistant  coatinqs.  '  Polymeric 
materials  (epoxies,  silicones,  polyimides  and  teflon)  with  thermal  capability  do  not  exhibit 
properties  of  tensile  strenqth,  elonqation  and  tear  strenqth  required  for  resistance  to  dynamic 
forces  of  impact.  Of  elastomeric  coatinqs,  polyurethane  on  aluminum  and  epoxy  laminate  sub¬ 
strates  was  found  to  be  up  to  5  times  as  resistant  to  subsonic  rain  impact  than  neoprene. 
Urethanes  in  which  the  polyol  molecular  weiqht  is  about  650  exhibit  best  resistance  In  simu¬ 
lated  rain  environment.  Polyurethane  coatinq  was  superior  when  exposed  to  rain  and  hall  in 
f I i qht  cond i t ions . 

Of  thermoplastics,  polycarbonate  polymer  (texan)  was  superior;  poorest  was  polytetra- 
fluoroethylene  (Teflon  TFE)  .  Of  metals,  the  stellite  alloys  (Co-Cr-W)  demonstrated  best 

(1,3) 

erosion  resistance.  Various  rainheads  of  stainless  steel,  neoprene,  solid  polycarbonate, 

polyurethane  and  metal  crossbar  designs  were  tested  for  PD  fuzes  at  various  sled  speeds  at 

Holloman.  For  example,  a  stainless  steel  sheet  15  mil  thick  was  punctured  after  2200  ft/sec 

(44) 

entry  into  the  2000  ft  lonq  Holloman  deluqe  rainfield.  The  forward  rim  surface  of  a  stain¬ 

less  steel  ogive  section  was  chewed  up  for  4600  ft/sec  entry.  Metals  are  not  marketly  strain 
rate  sensitive.  Impact  produces  plastic  deformation,  dependent  on  liquid  mass  qeometry,  when 
the  peak  pressure  exceeds  the  instantaneous  flow  stress  In  the  metal.  Polymers  and  glasses  are 
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strain  rate  sensitive,  and  tlie  tendency  is  to  elastic  deformation  in  the  areas  of  compression 
and  fracture  in  reqions  of  tension.  The  area  of  impact  is  found  relatively  undamaged,  tut 
surrounded  by  well  defined  rinq  cracks. 

(k5  k6) 

Water  which  flows  laterally  forms  a  protecting  liquid  film.  ’  An  increase  in 
the  density  of  drops  increases  the  water  film's  protective  role  and  reduces  deformation  from 
impact  shock.  If  the  film  is  too  thin,  multiple  reflection  in  the  thin  film  rapidly  brings  up 
pressure  to  full  value  for  "dry"  impact.  A  water  layer  comparable  to  half  a  cylinder  jet  ra¬ 
dius  acts  as  a  useful  cushion,  reducing  shock  pressure  by  factor  of  one  half. 

(A  7) 

In  an  investigation  of  the  splash  envelope  of  raindrops  at  low  impact  velocity 
(terminal  velocity  from  kO  ft  tower),  drops  qreater  than  about  2  or  3  mm  diameter  became  mea- 
sureably  flattened  on  the  bottom.  Drops  from  2.9  to  5-6  mm  fell  into  water  depths  of  0.1  to 
90  mm  above  a  smooth  plate  glass.  (A  splash  shape  was  not  formed  by  impact  on  smooth,  hard, 
dry,  horizontal  Surface;  rather  the  waterdrop  merely  spreads  horizontally  without  forminq  the 
splash  envelope  and  droplets  above  the  surface.)  The  splash  envelope  increased  to  its  greatest 
neight  when  the  water  depth  was  about  one-third  the  drop  diameter,  and  decreased  to  constant 
size  for  depths  greater  than  three  drop  diameters 

for  ceramics  and  high  hardness  metals.  It  is  thought  that  impact  fatigue  or  work 
Hardening  causes  small  imperfections  In  the  Surface,  to  be  removed  with  subsequent  impacts. 
Velocities  in  excess  of  Mach  1.3  (IkOO  ft/sec)  were  needed  to  damage  hign  density  alumina 
(AljOj)  with  2  mm  mercury  drops  However,  >t  was  reported  that  as  alumina  density 

increases  and  porosity  decreases,  the  bond  to  the  reinforced  laminate  becomes  weaker.  A  cor¬ 
rect  match  is  needed  for  the  thermal  expansions  of  laminate  and  coatlnq.  In  general,  the  alu¬ 
mina  coatings  were  about  2  1/2  times  more  rain  erosion  resistant  than  neoprene  coating.  It 
appears^^)  it  Is  not  possible  to  predict  rain  erosion  of  particular  ceramic  materials  at  high 
Mach  number  because;  (I)  theory  associated  with  prediction  of  forces  from  water  drop  impact 
on  solid  surfaces  is  insufficiently  developed;  (2)  property  data  for  ceramics  is  not  available; 
and  (3)  behavior  of  ceramics  under  multiple  impact  (fatigue)  is  not  understood.  Despite  the 


compressive  strenqth  of  alumina  ceramics,  they  were  damaqed  by  multiple  impacts  as  low  as 
820  ft/sec.  Pyroceram  9606  is  reported  superior  to  alumina  753- 
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Structures  (Figure  22) 


Investigations  indicated  that  spikes  ^  placed  on  missile  nosecones  effectively 

protect  radomes  from  erosion  at  velocities  up  to  3000  ft/sec  by  advancing  the  bow  shock  wave. 
To  protect  fragile  windows, an  open  end  frustrum  of  a  cone  was  attached  to  the  front  of  a 
missile  to  induce  raindrop  breakup.  The  captive  “air  cushion",  moving  at  essentially  flight 
speed,  provided  for  protection  up  to  Mach  2.5. 

FIGURE  22.  STRUCTURES  FOR  EROSION  ABATEMENT 

STRUCTURES:  (LOW  MACH  #) 

SP,KE  OPEN  FRUSTUM  (WALLEYE) 


Shock  Disintegration  pf  Props  (Mach  3  *  12) 


A  projectile  at  supersonic  velocity  through  air  with  a  blunt  nose  has  a  detached 
shock  wave  front;  a  sharp  point  projectile  has  a  shock  attached  to  the  point.  The  width  of 
the  shock  itself  is  exceedingly  small.  A  body  with  hemispherical  nosecone  of  radius  R  (ft) 
traveling  in  air  at  velocities  v  (ft/sec)  greater  than  the  speed  of  sound  (1090  ft/sec  at  0°F) 
is  preceded  by  a  detached  shock  wave  whose  stand  off  distance  d  (ft)  is  given  by:'**®^ 


d  -  R  (0. 129  ♦  7-61  x  l05/v2) 


SU'C* 

FUtMT 


(3M 


At  2000  ft/sec,  the  shock  wave  is  0.32,  and  at  6000  ft/sec  is  0.13  of  the  radius  away.  For  a 
body  of  revolution  with  a  flat  nose,  the  standoff  distance  is  about  0.8R  for  Mach  2,  and  about 
0.6R  for  Mach  6. 
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For  one  dimensional,  compressible,  Isentropic  perfect  gas  flow,  the  stagnation  tem¬ 
perature  Tj  (*R)  and  gas  density^  *n  *6e  *6ock  region  Is  given  by:^^ 

T 2  -  T,  (1  ♦  (k-l)  M2/2)  9  N  T,  (35) 

f2  "  Of, 

where  k  -  Ratio  of  heat  capacities  (Cp/Cv)  -  I.A 
M  -  Mach  number 

Pi ,  Tj  -  Atmospheric  density,  temperature  ("R) 

TABLE  21.  TEMPERATURE  (N)  AND  DENSITY  (D)  RATIO  FOR  NORMAL  SHOCKS 


Mach  Number  M 

N 

0 

1 

2 

1  . 

1.69 

1  . 

2.67 

3 

2.68 

3.86 

5 

5.8 

5.00 

10 

20.  A 

5.71 

The  reduction  of  the  drop  by  evaporation  would  be  a  function  of  rate  of  heat  transfer  and  time 
required  for  the  drop  to  traverse  the  shock  region.  This  mass  loss  by  heat  transfer  is  negli¬ 
gible  for  ordinary  rain  droplets,  but  the  cooling  effect  of  rain  helps  reduce  the  ablation 
rate  of  missile  material. 

A  projectile  or  missile,  dependino  on  its  design,  imparts  some  of  its  motion  to  the 
air  around  it.  For  hiqh  Mach  number,  up  to  a  million  G  acceleration  can  be  exDerienced  by 
drops,  and  if  the  Impact  with  the  front  air  stream  is  applied  far  enouqh  ahead  of  the  surface 
to  be  protected,  the  drop  may  be  shattered  into  droplets  loo  small  to  cause  damaqe.  When  the 
drop  is  exposed  to  an  airstream,  various  competing  forces  determine  whether  the  drop  will  be 
broken  up,  and  If  so,  the  time  for  the  process.  Forces  that  tear  apart  the  drop  scale 

approximately  with  dynamic  pressure  ^ VJ/7  and  the  tendency  to  maintain  drop  shape  with  iner¬ 
tia  forces,  surface  tension,  and  to  a  lesser  extent,  viscous  forces. 

A  criterion  to  predict  whether  or  not  a  liquid  drop  would  disintegrate  under  aero¬ 
dynamic  conditions  is  the  Weber  Number  W  (alternate  definition  has  7  in  denominator):  ratio  of 
inertial  energy  of  drop  to  its  surface  energy.  Some  find  the  critical  Weber  number  between 


W  ■  ^  u2  0/6 


(36) 


0  -  Drop  diameter 

S  -  Surface  tension  (.073  newtons/m  for  water  at  20° C ) 

2 ,  U2  -  Air  density,  speed  relative  to  drop 

Quantitative  results  of  raindrop  breakup  effects  (Fiqure  23)  from  a  strong  shock 

FIGURE  23-  SHOCK-UATERDROP  COLLISION  IN  LAB  COORDINATES 


(hach  3  to  12)  were  obtained  in  shock  tubes  with  drop  diameters  from  0.5  to  2.5  nm  with 


pressures  from  IA0  to  760  torr.  The  displacement  *  after  actual  time  t  of  shock  front  passage 
of  an  initially  motionless  drop  followinq  contact  with  a  shock  front  is  reported  as: 

«  ~0. 8  0  T*  (?7) 


diameter  T  =  tU2\[p 
time  Q  N 


where  D  -  Oriainal  drop  di 
T  -  Dimensionless 
^2'  Pl  "  Gas  density  in  shock  repion,  liquid  density 

U2  -  Flow  velocity  in  shock  reqion  in  lab  coordinates 
x  -  Distance  traveled  by  drop  In  lab  coordinates 

Equation  37  corresponds  to  a  constant  acceleration  path  with  an  averaqe  draq  coef¬ 
ficient  of  about  2.1  based  on  the  oriqinal  drop  cross  section.  Table  22  lists  standard  (NACA 
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Report  1235)  pressure-duns  I ty-temperature  values  for  various  altitudes. 


TA81E  22.  PRESSURE -DENS I TY-TEMPERATURE  FOR  VARIOUS  ALTITUDES 


Alt! tude 
(ft) 

Atmospheres 

Pressure 

(torr) 

Pressure 

(lb/fn2) 

Density 

(grams/cc)  x  10’ 

Temperature 
(*  K) 

0 

1. 

760 

Ik. 69 

12.3 

288.2 

10,000 

.687 

523 

10. 1 

9.07 

268.3 

20,000 

.  k6 

350 

6.75 

6.5k 

2k8 .6 

30,000 

.297 

226 

k.  36 

k.58 

228.7 

kO ,000 

.185 

Ik  1 

2.72 

3.01 

216.6 

50,000 

.Ilk 

87 

1  .68 

1.86 

216.1 

60,000 

.071 

5k 

1 .0k 

1.15 

216.6 

100,000 

.011 

6 

0.158 

0.17k 

216.6 

The  drop  is  always  flattened  Iritlally  by  the  pressure  differential  from  air  flow. 


There  are  four  descriptive  modes  for  the  disruption  of  the  liquid  drop  by  alrstream  Impact: 

(I)  Drop  oscillates  until  division  in  tv/o,  at  low  speeds;  (2)  At  sliqhtly  higher  speeds,  the 
drop  is  severely  distorted  into  parachuti  shape  and  soon  shatters;  (3)  At  still  higher  speed, 
continuous  stripping  mode  of  breakup  occurs  In  which  a  spray  is  formed  at  the  periphery  of  the 
drop  and  swept  Into  the  wake;  (k)  At  extreme  hiqh  speeds  (catastrophic)  ,  drops  shatter  very 
rapidly  in  a  distinct  mode. 

The  transition  from  stripping  to  catastrophic  mode  occurs  early  in  the  shock  tube 
tests  at  Mach  II,  at  intermediate  time  at  Mach  6,  and  at  later  time  (if  at  all)  at  Mach  3. 

X-ray  data  suggested  time  for  complete  breakup  due  to  stripping  has  a  constant  value  T  &  3.5 
or  *  s  10  D,  independent  of  Mach  number.  The  mass  m  of  the  drop  remaining  of  the  initial  mass 
ny,  as  a  function  of  time  could  be  correlated  roughly  bv  the  formula: 

m  -  m<j  (1  ♦  cos  (TIT/3. 5)  )/2  (39) 

Uhen  the  critical  Weber  number  Is  exceeded,  the  first  mode  occurs.  For  I  to  2  mm  dia¬ 
meter  drops,  the  critical  value  ranges  from  about  k  to  13,  corresponding  to  velocities  of  kO  to 
200  ft/sec  at  sea  level .  ^  As  W  Is  Increased,  there  is  a  small  range  of  values  for  which  the 
second  mode  occurs.  For  large  values  (about  3  times  the  critical),  breakup  occurs  by  flattening 
and  stripping.  For  high  velocity  flight,  the  third  mode  Is  relevant  wherein  the  creation  of 
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small  fraqments  leads  to  Increase  of  surface  to  volume  ratio,  and  sudden  increase  in  accelera¬ 
tion.  The  dimensionless  time  T  elapse  between  shock  contact  and  drop  breakup  in  the  catastro¬ 
phic  mode  is  correlated  by  equation  40: 

-0.25 

T  -  45  W  (40) 

An  i nvest ioat Ion  of  a  2.7  mm  diameter  water  drop  at  Mach  1.5  indicated  substantial 
deformation  several  inches  beyond  the  shock  front  passaqe.^'^ 

Waterdrops  will  be  shattered  as  they  pass  throuqh  the  bow  shock  region  that  precedes 
a  body  travelling  at  supersonic  velocity  only  if  the  time  of  traversal  in  the  shock  region  is 
greater  than  the  break-up  time.  The  traverse  time  is  determined  by  the  velocity  and  geometry 
of  the  leading  surface  of  the  body.  Data  indicates  that  raindrops  of  normal  size  will  impact 
unshattered  on  the  hemispherical  nosecone  of  a  missile  at  1000  ft/sec  if  the  nosecone  radius 
is  less  than  about  0.7  ft;  at  2000  ft/sec  if  the  radius  is  less  than  about  1.6  ft;  and  at 
10,000  ft/sec  if  the  radius  is  less  than  2.8  ft.^^  It  is  expected  that  raindrops  will  not 
be  shattered  in  passing  through  the  bow  shock  layer  of  most  missiles,  and  serious  erosion 
would  occur  with  most  nosecone  materials  on  missiles  havinq  nosecone  radii  less  than  1.5  ft  at 
velocities  qreater  than  about  3000  ft/sec,  with  qreater  damaqe  caused  by  those  drops  strikinq 
normal  to  the  surface. 


Hypervelocity  Impact 

Rather  than  erosion,  craterinq  from  extreme  compressive  force  of  drops  becomes  the 
dominant  form  of  damage  at  velocities  above  a  few  Mach  numbers.  As  the  velocity  approaches 
Mach  10  to  20,  water  drops  cause  the  same  penetration  damage  as  a  ductile  solid  of  the  same 
density,  as  the  only  appreciable  parameters  are  velocity  and  density.  A  jet  of  water  at  5000 
ft/sec  has  sufficient  force  to  cut  metal.  Data  for  estimating  drop  damaqe  at  velocities 
qreater  than  about  3000  ft/sec  is  obtained  from  impacts  of  single  drops  of  liquids  or  solid 
particles  on  semi  -  Inf i n 1 te  tarqets  (with  thickness  at  least  five  times  penetration  depth). 
Damaqe  is  normally  reported  as  depth  of  penetration  (P)  divided  by  drop  diameter  ( D) .  Figure 
24  shows  the  crater  profile,  though  the  crater  top  has  a  jagged,  peeled  over  edge. 
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FIGURE  2 It.  CRATERING  CHARACTERIZATION 


For  velocities  from  5000  to  35,000  ft/sec,  the  followino  equatinn^-*^  Is  proposed  to  explain 
87  percent  of  total  variance  of  data  in  this  ranoe  : 

n  .673  .449 

0-33  t_P _ V 

p  -  0.772  (V  )  .1*26  .275  (42) 

p  S 


Another  equation 


(58) 


for  hypervelocl ty  impact  that  correlates  a  wide  variety  of  metal -to-metal 


impact : 


P  -  0.31 1  0 

>p' 

.167 

1  >i 
0.1 

.  s  J 

.282 


(43) 


Cal culat ions ^ ^  with  HELP  hydro  code  (mul t i -mater ia I  Eulerian,  compressible  fluid, 
elastic-plastic  flow,  time  and  20  space)  with  a  3-4  mm  water  drop  onto  606I-T65I  aluminum  tar- 
qet  indicated  collision  lasts  approximately  5  microseconds,  with  some  tarqet  rebound  for  the 
4655  ft/sec  encounter.  Results  are  in  Table  23  where  Q  is  the  distance  P  plus  heiqht  of  the 
deformation  above  the  surface  line,  with  some  extrapolated  results  from  a  water-qelatin  pro¬ 
ject  i I e  experiment . 


TABLE  23.  HELP  2D  C00E  CALCULATIONS  AND  EXTRAPOLATED  WATER-DROP  (3.4  run)  EXPERIMENT 


|  Crater  Dimensions  (Calculated) 

Extrapolated  Experiment 

( Inch) 

(Inch) 

Velocity  (ft/sec) 

P 

Q 

w 

P 

W 

4655 

.018 

.048 

.194 

.012 

.18 

9180 

.072 

.126 

.278 

Water  produces  a  very  shallow  and  flat-bottomed  crater  at  5000  ft/sec  but  the  crater  is  nearly 
hemispherical  at  10,000  ft/sec  and  has  nearly  the  same  depth  as  for  an  aluminum  impact  of  the 
same  enerqy. 


Miscel lanea 

(59) 

Initial  studies'  showed  similarities  In  the  behavior  of  hail  to  that  of  water 
drops  durinq  the  Impact  phase.  Hall,  spherical  and  bullet  shaped,  in  diameters  from  12.7  to 
25-4  mm  were  fired  from  a  hlqh  pressure  qas  qun  at  speeds  up  to  3000  ft/sec  onto  aluminum  alloy 
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(0.6  to  3  mm)  and  polymetnylmethactylate  (6 .35  to  <  j 
Photoqraphs  showed  that  shortly  after  impact,  a  t  .vi  -  ■  '  •  i  .  :s  tmii..  a. 

there  were  no  sions  of  spaMinq  from  the  rear  Sir  ace  pf  roicr-ii-  *  ‘  qro  ' 

taiven  as  0.9.  The  shock  wave  velocity  in  ice  is  about  twice  oat  in  water  iublt.:,  i'i 
list  reported^^)  hail  and  ice  crystal  distributions.  Hail  can  exist  at  altitudes  up  t  V' 
60,000  ft. 


TABLE  26.  POSSIBLE  HAILS  fONt  SIZE  Jl  STI' I  B6T  I  OH 


Maximum  Hailstone  Diameter  (Centimete: 


Diameter  R  art  ae  (cm) 

}  cm  6  cm 

No.  Ha i 1 stones/Cuh i c 

1  0  cm 

He  ter 

0-1 

22.* 

9  ■  * 

/.* 

1-2 

3.8 

0.35 

0  27 

2-3 

3.17 

0.075 

0  356 

3-6 

0.027 

0.021 

6-5 

— 

0.01 3 

0  .  u  1 C 

5-6 

0.0007 

C  0-.  5 

6- 7 

7- 8 

0.003 

0.002 

a  vj 

... 

0.0CI5 

9-10 

0.00! 

itM.iivaient  ft(mr*/hr) 

>00 

500 

!  (mm6/m3)  10  cr  radat 

S  J  X  I  L." 

5.6  x  10 

I Equi valent  L  (o/m3) 

3.8 

33 

6.3 

1*  Probably  a  considerable  nnderes t :  ■■  1 1 « 

f  f  the  sma 1 1 

: s  t  stones  in 

this  ranqe. 

TABLF  25.  ICE  (b, STAI  CONTENT  OF  CLOUD 


Air  Temperature 
CC) 

Altitude  Ranoe  (ft) 

Ice  Crystal  Content 
(Oram/m3) 

Horizontal  Extent 
(Miles) 

0  to  -20 

10,000  -  30,000 

8.0 

0.5 

(3  -  9  km) 

5.0 

3.0 

2.0 

50. 

1.0 

300. 

-20  to  -60 

15,000  -  60,000 

5.0 

3. 

(6.5  -  12  km) 

2.0 

10. 

1  .0 

50. 

_ _ _  . 

0.5 

300. 

Over  temperature  ranqe  0  to  -I0*C,  assume  that  ice  crystals  may  be  mixed  with  water 
drops  (with  maximum  diameter  of  2  ntn)  to  up  to  content  of  1  qm/m3  or  half  the  total  content, 
whichever  is  less;  the  total  content  remainlnq  unehanned  Mean  diameter  mav  be  taken  as  I  mm. 


Precipitation  has  been  measured  by  use  of  radar  reflectivity  factor  2  which  Is  di¬ 
rectly  proportional  to  ^NjDj6  which  in  turn  is  related  to  rainfall  rate  R.  For  example,  one 
relation^  for  all  storms  Is  Z(mm6/m3)  »  J72  R(mri/hr)  *  ,J,7. 

A  non-metal  lie  surface  can  acquire  a  sliqhtly  negative  charge  from  rain  impact. 

Values  of  charge  per  drop  have  been  cited^'^  as  0.37  (cloudburst)  to  0.0069  (moderate  rain) 
esu/drop. 

(A  copy  of  the  papers  presented  at  the  Third  International  Conference  on  Rain  Erosion 
and  Associated  Phenomena  (1970)  was  not  available  at  the  time  of  this  writing.) 
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